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Preface 


The  purpose  of  this  study  was  to  develop  a  control  system  capable  of  controlling  a 
formation  of  aircraft.  Special  attention  was  placed  on  the  analytical  analysis  of  the  control 
problem  and  the  development  of  accurate,  yet  simple  models  of  the  formation  aircraft. 
Motivated  to  improve  the  mission  effectivness  of  the  Air  Force  Special  Operations  Forces, 
the  C-130H  aircraft  was  chosen  as  the  primary  formation  aircraft. 

The  analytical  analysis  was  performed  on  the  formation  flight  control  problem  clearly 
indicates  the  necessity  of  integral  control  in  the  design  of  the  automated  formation  flight 
control  system.  This  analysis  also  revealed  the  interelationship  of  the  formation  geometry, 
aircraft  models  and  controller  parameters  on  system  stability.  Further  investigation  led  to 
the  design  of  a  three  dimensional  formation  control  system,  two  wing  aircraft  simulations 
and  the  development  of  a  visualization  tool  in  evaluating  the  formation  flight  control 
system. 

In  performing  the  research  and  writing  this  thesis,  I  have  had  a  great  deal  of  help 
from  others.  I  am  deeply  indepted  to  my  thesis  advisors.  Dr.  Meir  Pachter,  Dr.  John 
J.  D’Azzo  and  Dr.  C.  Houpis.  Their  continued  guidance  and  technical  assistance  proved 
invaluable  during  difficult  times.  Specifically,  I  am  grateful  for  their  mathematical  and 
literary  expertise.  I  am  also  grateful  for  the  long-distance  support  and  encouragement 
given  by  my  parents  throughout  this  endeavor.  Their  concern  and  interest  was  extremely 
beneficial. 


Louis  E.  Buzogany 
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Abstract 

In  this  study,  an  automated  formation  control  system  for  an  aircraft  formation  com¬ 
prised  of  one  lead  and  multiple  wing  aircraft  is  analyzed.  Second-order  models  of  the  C-130 
aircraft  are  developed  in  order  to  accurately  model  the  flying  qualities  of  large  aircraft. 
This  automated  formation  control  system  is  capable  of  controlling  the  C-130  aircraft  in 
maneuvering  formation  flight,  thus  reducing  the  wing’s  pilot  workload. 

During  formation  flight,  the  wing  aircraft  continuously  measures  the  lead  aircraft’s 
relative  position  with  an  ideal  on-board  position  sensor.  This  information,  in  addition  to 
Proportional  Plus  Integral  feedback  control,  is  used  to  maintain  the  aircraft  in  formation. 
The  control  of  each  wing  aircraft  is  assumed  to  be  independent  of  other  wing  aircraft. 
Other  than  for  nominal  formation  separation  commands,  no  continuous  communication  is 
assumed  between  the  formation  aircraft. 

An  analytical  analysis  of  the  formation  control  problem  reveals  that  integral  control 
is  needed  to  achieve  zero  steady  state  error  in  the  separation  distances  (after  a  forma¬ 
tion  maneuver  is  executed).  This  conclusion  is  confirmed  using  computer  simulations.  An 
analytical  method  of  selecting  the  Proportional  Plus  Integral  parameters  is  developed  by 
identifying  the  dominant  system  dynamics  and  residues  of  the  step  response.  In  an  attempt 
to  reduce  the  fuel  consumption  of  the  wing  surcraft  during  formation  heading  change  ma¬ 
neuvers,  an  alternate  control  system  is  designed  to  conserve  the  energy  of  the  wing  aircraft. 
Thus,  an  altitude  change  is  used  to  control  the  wing  aircraft’s  velocity,  minimizing  the  use 
of  the  throttles  during  a  heading  change  maneuver. 

The  resulting  automated  formation  control  system  effectively  maintains  the  formation 
of  aircraft  through  a  combination  of  velocity,  heading,  and  altitude  changes.  There  is  zero 
steady  state  error  for  all  maneuver  and  separation  distance  change  responses.  Adflit.intia.lly, 
the  separation  transients  axe  such  that  no  collisions  occur  between  the  aircraft.  Both 
analytical  equations  and  computer  simulation  generated  time  responses  are  provided  for 
reference. 
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AUTOMATED  CONTROL  OF  AIRCRAFT  IN  FORMATION  FLIGHT 


/.  Introduction 


1.1  Overview  of  the  Thesis 

The  introduction  and  general  background  information  relevant  to  the  formation  flight 
control  problem  is  provided  in  Chapter  I.  The  simulation  development,  aircraft  autopilot 
models,  and  formation  controller  design  is  presented  in  Chapter  II.  A  method  of  evaluating 
the  transient  response  of  the  formation  flight  control  system  by  displaying  the  lead  air¬ 
craft’s  trajectory  in  the  wing  aircraft’s  reference  frame  is  presented  in  Chapter  III.  A  two 
wing  aircraft  formation  and  performance  simulation  is  presented  in  Chapter  IV.  A  detailed 
derivation  of  a  linear  and  parametrized  Multi-Input  Multi-Output  (MIMO)  plant  used  to 
analytically  determine  the  optimal  Proportional  Plus  Integral  (PI)  control  parameters  for 
the  formation  flight  control  problem  is  presented  in  Chapter  V.  Chapter  VI  presents  the 
theoretical  analysis  and  computer  simulation  of  an  3-Dimensional  (3D)  formation  flight 
control  system  based  on  energy  conservation  of  the  wing  aircraft.  A  performance  eval¬ 
uation  of  the  formation  flight  control  system  incorporating  second-order  aircraft  models 
is  presented  in  Chapter  VII.  The  method  used  to  obtain  second-order  aircraft  models  is 
presented  in  Appendix  A.  The  design  and  computer  simulation  of  longitudinal  and  lateral 
autopilots  developed  around  the  C-130  aircraft  is  included  in  Appendix  B.  Conclusions 
and  areas  requiring  further  study  are  presented  in  Chapter  VIII. 

1.2  Background 

Military  pilots  are  required  to  fly  many  missions  that  require  extended  periods  of 
formation  flight.  During  these  periods,  the  pilot’s  attention  is  divided  between  monitoring 
the  instruments  of  his  aircraft,  and  maintaining  the  correct  separation  distance  to  the 
adjacent  aircraft.  The  denser  the  formation  of  aircraft,  the  more  time  the  pilot  spends  on 
maintaining  the  formation.  As  other  conditions,  such  as  visibility,  navigation  and  enemy 
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threat  begin  to  require  more  of  the  pilot’s  attention,  the  more  difficult  it  is  for  him  to 
maintain  the  formation(2)[l'2]. 

An  automated  formation  ffight  control  system  is  developed  to  relieve  the  strain  on 
military  pilots  during  formation  ffight.  This  system  is  analogous  to  an  autopilot  used  on 
a  single  aircraft,  except  that  this  "autopilot”  is  used  to  control  a  formation  of  aircraft. 

The  Air  Force  Special  Operations  Forces  (SOF)  perform  missions  involving  extended 
periods  of  formation  ffight  that  place  heavy  workloads  on  the  pilots.  SOF  pilots  are  tasked 
to  conduct  overt,  clandestine,  or  covert  missions  which  can  range  from  routine  training 
missions  to  highly  sensitive  missions  of  national  importance(2)[l-2].  Many  times  these 
missions  require  close  formation  ffight  at  low  altitudes,  using  similar  or  dissimilar  aircraft. 
The  C-130  mrcraft  are  predominantly  used  during  these  missions. 

The  C-130  aircraft,  models  C-130A  and  C-130B,  fly  several  typical  formations,  in¬ 
cluding  the  trail  and  diaunond  formations.  The  trail  formation  is  shown  in  Fiqure  1.1, 
while  the  diamond  formation  is  shown  in  Fiqure  1.2.  The  trail  formation  is  important 
since  a  minimum  amount  of  lamd  maiss  is  overflown  by  the  formation  aircraft,  translating 
to  a  reduced  probability  of  detection  by  ground  forces.  The  diamond  formation  maintains 
good  visibility  between  the  various  aircraft. 
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Various  maneuvers  are  performed  while  the  aircraft  are  flying  in  formation.  These 
maneuvers  are  initiated  by  the  lead  aircraft,  without  prior  knowledge  of  the  wing  aircraft. 
The  lead  aircraft  may  initiate  a  change  in  velocity,  heading,  altitude,  or  a  combination  of 
the  three.  Fiqure  1.3  shows  a  heading  change  maneuver  and  Fiqure  1.4  shows  an  altitude 
change  maneuver.  These  maneuver’s  are  executed  from  a  diamond  formation  and  a  tra.il 
formation  respectively.  Depending  on  varying  conditions  and  mission  requirements,  the 
formation  used  by  SOF  aircraft  may  change  throughout  the  flight.  A  change  from  trail 
formation  to  diamond  formation  is  shown  in  Figure  1.5. 
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Figure  1.3.  Diamond  Formation  Heading  Change  Maneuver 
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Figure  1.4.  Trail  Formation  Altitude  Change  Maneuver 
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Figure  1.5.  Trail  to  Diamond  Formation  Change 

The  most  difficult  maneuver  for  the  wing  aircraft  is  a  large  heading  change  in  a 
diamond  formation.  This  type  of  maneuver  is  shown  in  Figure  1.3.  The  wing  mrcraft 
must  initially  increase  its  velocity  and  traverse  the  outside  path  in  order  to  track  the  lead 
aircraft,  following  which  the  wing  aircraft  will  reduce  its  speed.  This  modulation  of  the 
velocity  results  in  increased  fuel  consumption.  In  an  attempt  to  save  fuel,  it  is  customary 
for  the  pilot  of  the  wing  aircraft  to  initially  descend  to  increase  his  velocity  during  the 
maneuver,  and  then  to  ascend  to  bleed  off  any  excess  airspeed.  Considering  a  single  wing 
aircraft  formation,  the  wing’s  pilot  may  also  initiate  an  altitude  change  in  addition  to  a 
heading  change.  An  example  of  such  a  maneuver  is  shown  in  Figure  1.6.  This  type  of 
maneuver  is  addressed  in  detail  in  Chapter  VI. 

1.3  Problem  Statement 

The  problem  of  improving  SOF  mission  effectiveness  through  the  implementation  of 
an  automated  formation  flight  control  system  is  addressed  in  this  thesis.  The  automated 
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Figure  1.6.  Diamond  Formation,  90  Degree  Heading  Change  using  Energy  Conservation 
and  Formation  Transpositon 

formation  flight  control  system  will  both  relieve  pilot  fatigue  and  allow  the  performace  of 
other  tasks  during  the  formation  flight.  Accurate  aircraft  models  and  thorough  control 
system  analysis  directly  impacts  the  validity  of  the  final  results.  Thus,  the  purpose  of  this 
thesis  is  to  develop  models  and  methods  needed  to  design  an  automated  formation  flight 
control  system  to  control  a  formation  of  aircraft. 


1.4  Summary  of  Current  Knowledge 

Two  previous  theses  research  efforts  by  Capt  Paul  Rohs  and  Capt  John  Dargan  at  the 
Air  Force  Institute  of  Technology  (AFIT)  investigated  the  development  of  an  automated 
formation  flight  control  system.  Both  used  a  digital  computer  simulation  for  development 
and  evaluation.  Capt  Rohs  developed  a  control  system  capable  of  controlling  a  formation 
of  similar  or  dissimilar  aircraft  through  three  separate  maneuvers(9).  These  maneuvers 
included  a  heading  change,  a  terrain  avoidance  altitude  change,  and  a  formation  change. 
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Aircraft  position  was  measured  using  a  Cartesian  reference  system.  It  was  assumed  that 
perfect  lead  aircraft  position  information  was  available  to  each  wing  aircraft  control  system. 
First  order  mathematical  models  of  different  aircraft  were  used  in  the  simulations.  Capt 
Dargan  continued  Capt  Roh’s  research  and  developed  a  control  system  composed  of  a  PI 
controller  and  a  linear  mixer.  Formations  consisting  of  one  lead  and  one  wing  aircraft  were 
used  for  simulation  and  evaluation.  First  order  models  of  C-130  aircraft  were  also  used. 
The  results  of  this  research  showed  that  a  PI  controller  can  be  developed  for  an  automated 
formation  flight  control  system  that  maintains  formation,  without  collisions,  and  provides 
satisfactory  transient  behavior  with  zero  steady  state  error(2)[6-5]. 

1.5  Research  Objective 

The  objective  of  this  research  is  to  continue  the  development  of  an  automated  for¬ 
mation  flight  control  system  through  an  analytical  analysis  of  the  formation  flight  control 
problem  and  computer  simulations  using  more  realistic  aircraft  models. 

1.6  Research  Questions 

The  following  questions  are  answered  in  this  thesis. 

1.  Can  second-order  aircraft  models  be  developed  in  order  to  more  accurately  model 
the  flight  performance  of  the  C-130  aircraft? 

2.  What  is  the  effect  of  second-order  aircraft  models  on  the  response  of  the  formation 
flight  control  system. 

3.  Can  the  current  formation  control  system  containing  the  PI  controller  and  linear 
mixer  be  used  to  address  the  two  wing  aircraft  formation  flight  control  problem? 

4.  Can  the  derivation  of  a  linear  and  parametrized  MIMO  plant  be  used  to  analytically 
determine  the  optimal  control  gains  for  the  formation  flight  control  problem.? 

5.  Can  a  flight  control  system  be  developed  to  perform  3-Dimensional  maneuvers  in 
order  to  reduce  the  fuel  consumption  of  the  wing  aircraft  during  heading  change 
maneuvers? 
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6.  Can  a  tool  for  evaluating  the  formation  flight  control  system  be  developed  by  graph¬ 
ically  displaying  the  trajectory  of  the  lead  aircraft  with  respect  to  the  wing? 

1.7  Assumptions 

In  support  of  the  thesis  objectives,  and  questions  listed  earlier,  the  following  assump¬ 
tions  are  made: 

The  lead  aircraft  is  piloted  by  a  human  pilot,  while  the  wing  aircraft  are  controlled  by 
the  formation  control  system.  The  formation  of  aircraft  is  referenced  to  the  lead  aircraft. 
Thus  all  formation  maneuvers  are  initiated  by  the  pilot  of  the  lead  aircraft  with  no  prior 
knowledge  of  the  maneuvers  given  to  the  wing  aircraft.  Each  wing  aircraft  contains  the 
formation  flight  control  system  avionics,  independent  of  the  system  mounted  on  the  other 
aircraft  in  the  formation. 

Each  wing  aircraft  is  assumed  to  possess  an  on-board  sensor  capable  of  providing 
precise  position  information  relative  to  the  wing  aircraft 12].  This  sensor  measurement 
data  is  used  to  track  the  lead  aircraft. 

Each  aircraft  in  the  formation  contains  an  automated  flight  control  system  capable  of 
maintaining  or  changing  the  aircraft’s  velocity,  heading  and  altitude.  Thus,  each  aircraft  is 
endowed  with  an  Mach- Hold  autopilot,  capable  of  adjusting  the  aircraft’s  forward  velocity 
to  control  separation  along  the  flight  path,  a  Heading-Hold  autopilot  for  performing  coor¬ 
dinated  turns  for  lateral  separation  control,  and  an  Altitude-Hold  autopilot  for  climbing  or 
ascending  for  vertical  separation  control.  The  mathematical  autopilot  models  used  in  this 
simulation  are  based  on  the  performance  of  the  C-130  aircraft.  Both  first  and  second  order 
aircraft/autopilot  models  are  used.  The  first  order  models  are  used  in  analytical  analyses 
and  for  performance  comparison  against  previous  research. 

The  initial  conditions  for  all  formation  maneuvers  are  straight  and  level  flight  in  a 
constant  formation.  Formation  maneuvers  are  executed  one  at  a  time;  however,  simulta¬ 
neous  heading  change  maneuvers  and  formation  transpositions  arc  considered. 


1-8 


1.8  Scope 

The  objective  of  this  research  is  to  develop  a  mathematical  control  law  capable  of 
controlling  a  formation  of  aircraft  under  realistic  (simulation)  conditions.  The  system 
designed  for  this  research  is  developed  through  the  use  of  the  MATRIXx  CAD  tool.  Even 
though  two  different  models  of  the  C-130  are  available,  the  C-130B  is  used  as  the  primary 
formation  aircraft.  Under  the  assumption  that  perfect  position  information  is  available, 
position  sensor  models  are  not  used.  This  research  is  not  intended  to  result  in  the  hardware 
fabrication  of  the  formation  control  system.  Hence,  hardware  specifications  such  as  power, 
weight  and  size  requirements  are  not  addressed(2)[l-13]. 

1.9  Standards 

Conventional  control  system  analysis  and  design  standards  are  used  in  this  research. 
Control  system  responses  to  input  commands  are  analyzed  with  respect  to  stability,  tran¬ 
sient  behavior  and  steady  state  values.  Time-response  specifications  of  rise  time,  settling 
time,  peak  overshoot  and  steady  state  error  is  emphasized  (7:90).  While  simulation  tools 
are  available  that  are  capable  of  measuring  these  time  response  specifications,  the  primary 
means  of  evaluating  the  transient  response  is  the  shape  of  the  output  waveform. 

1.10  Approach/Methodology 

The  formation  control  systems  developed  by  Captains  Rohs  and  Dargan  serve  as 
q  foundation  for  this  research.  The  PI  controller  developed  by  Capt  Dargan  is  verified 
through  simulation  on  the  digital  computer.  Upon  verification,  this  system  is  used  as  a 
starting  point,  and  modified  as  needed  during  the  design  process. 

The  premise  of  using  Dargan ’s  formation  control  system  as  a  starting  point  is  justified 
since  his  research  showed  that  the  system  worked  under  the  given  assumptions.  Therefore, 
since  a  system  has  been  previously  designed  and  has  been  shown  to  work,  it  is  not  necessary 
to  start  over  and  design  a  completely  new  system.  This  system  is  used  for  the  majority  of 
the  simulations  and  theoretical  analyses. 
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1.11  Benefits  of  the  Research 


This  research  provides  a  method  of  designing  a  control  system  able  to  maintain  a 
formation  during  maneuvers  and  able  to  command  separation  distances.  Related  research 
is  being  conducted  by  the  Intra  Formation  Positioning  System  (IFPS)  program  at  the  Flight 
Dynamics  Laboratory.  This  thesis  provides  an  additional  approach  to  the  formation  flight 
control  problem.  The  ideas  and  information  from  this  research  may  be  used  in  the  IFPS 
program  to  help  solve  a  real  world  problem. 

1.12  Literature  Review 

The  following  technical  literature  dealing  with  the  application  of  control  theory  to 
aircraft  formation  flight  is  reviewed.  This  review  includes  the  equation  of  Coriolis,  Porter’s 
Method  for  controller  design,  and  large  aircraft  flying  qualities.  This  review,  in  addition 
to  the  previous  research  by  Dargan  and  Rohs,  gives  the  background  information  needed  to 
develop  the  simulation  model  of  a  formation  control  system. 

1.12.1  Equation  of  Coriolis  The  equation  of  Coriolis  is  needed  when  the  motion 
of  one  object,  relative  to  another  object  in  a  different  rotating  reference  frame,  is  desired. 
The  motion  of  this  object,  as  viewed  from  the  rotating  reference  frame,  consists  of  the 
motion  as  seen  from  its  moving  frame,  plus  the  motion  resulting  from  the  relative  angular 
velocity  of  the  moving  frame  with  respect  to  the  inertial  reference  frame(l)[489].  This 
situation  is  illustrated  in  Figure  1.7.  The  equation  of  Coriolis  is 

R,- =  Rp+a;,  p  xRp  (1.1) 

where 


Rj  =  velocity  of  object  in  frame  i 

Rp  =  velocity  of  object  in  frame  p 

wip  =  angular  velocity  of  frame  p  with  respect  to  frame  i 
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Figure  1.7.  Inertial  and  Rotating  Frames  of  Reference 


Rp  =  position  of  object  in  frame  p 

The  equation  of  Coriolis  is  also  referred  to  as  the  velocity  transformation  law.  This  equation 
is  useful  in  expressing  the  lead  aircraft’s  velocity  in  the  wing  aircraft’s  rotating  reference 
frame. 

1.12.2  Porter’s  Design  Method  The  decoupling  of  the  outputs  is  an  essential 
requirement  in  many  multi-input  multi-output  (MIMO)  control  systems.  Some  design 
methods  require  that  all  the  states  must  be  fed  back  in  order  to  achieve  output  decoupling. 
Depending  upon  the  dimension  of  the  state  vector,  it  may  not  be  physically  possible  to 
include  enough  sensors  which  will  measure  all  of  the  states.  It  may  be  necessary  in  this 
case  to  use  a  state  estimator  or  observer  to  reconstruct  the  states,  and  then  feed  back  these 
estimated  states.  However,  a  control  design  method  which  uses  only  output  feedback  to 
generate  an  error  vector  avoids  the  requirement  for  measuring  or  reconstructing  the  entire 
state  vector.  One  such  method  is  that  developed  by  Brian  Porter  (7)  (8).  His  method 
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of  designing  a  high  gain  proportional  plus  integral  controller  produces  output  decoupling 
and  leads  to  very  fast  tracking  of  the  command  input  by  the  output.  The  MIMO  plant  is 
represented  by  the  standard  state  and  output  equations,  of  the  respective  forms: 


x(t) 

=  Ax(t)  +  Bu(t) 

(1.2) 

y(t) 

=  Cx(t) 

(1.3) 

where 

t  =  time 

x(t)  =  state  vector 
A  =  state  transition  matrix 
B  =  input  matrix 
u(t)  =  input  vector 
y{t)  =  output  vector 
C  =  output  matrix 

A  requirement  for  Porter’s  design  method  is  that  the  matrix  product  CB  have  full  rank. 
When  CB  has  fuU  rank  the  plant  is  described  as  “regular”.  When  CB  does  not  have 
full  rank  the  plant  is  described  as  “irregular”.  In  the  case  of  regular  plants  the  controller 
implements  a  proportional  plus  integral  control  law  in  the  forward  path  of  the  control 
system.  In  the  case  of  irregular  plants,  the  proportional  plus  integral  control  is  augmented 
with  an  inner-loop  which  provides  extra  measurements  for  control  purposes  (3:660-661). 

1.12.3  Large  Aircraft  Flying  Qualities  Flying  qualities  are  defined  as  those 
characteristics  of  an  aircraft  that  govern  the  ease  and  precision  with  which  a  pilot  is  able 
to  perform  his  mission  (5)[456].  There  are  significant  differences  in  the  flying  capabilities 
between  large  and  relatively  small  aircraft.  These  differences  are  reflected  in  the  short 
period  frequency  in  the  longitudinal  mode,  the  time  to  bank  in  the  lateral  mode,  and  the 
time  delay  for  all  axes.  Of  special  importance,  in  terms  of  modelling  the  flying  performance 
of  the  aircraft,  is  the  time  delay. 
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The  time  delay  is  measured  as  the  time  from  the  initiation  of  a  step  control  input, 
until  the  first  indication  of  the  overall  aircraft  response  (5)[459].  Even  though  time  de¬ 
lays  degrade  the  flying  qualities  rating  of  the  aircraft,  and  may  result  in  Pilot  Induced 
Oscillations  (PIOs),  time  delays  greater  than  1  second  have  been  measured  (5)[460].  The 
presence  of  these  delays,  especially  in  large  aircraft,  reinforces  the  motivation  to  develop 
second  order  models  for  the  formation  control  system  simulation. 

1.13  Materials  and  Equipment 

The  following  equipment  and  materials  are  required  for  this  research: 

•  The  MATRIXx  Computer  Aided  Design  tool.  Both  the  Unix  Open  Windows  and 
the  IBM  PC  compatible  versions  are  used. 

•  The  Sun  SPARC  workstation  computers  located  in  the  AFIT/ENG  computer  lab  in 
room  133  in  building  640. 

•  Both  personal  and  AFIT/ENG  IBM  compatible  personal  computers. 

•  Generic  CADD  5.0  computer-aided  drafting  software,  personal  copy. 

•  Latex  document  preparation  software,  Unix  and  IBM  PC  versions. 
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II.  Simulation  Development 


A  block  diagram  representing  the  computer  simulation  used  to  model  the  formation 
flight  control  system  is  shown  in  Figure  2.1.  The  lead  aircraft  is  totally  independent  of 
the  wing  aircraft.  The  information  about  the  lead  aircraft  available  to  the  wing  aircraft  is 
acquired  by  its  external  sensors.  The  external  sensors,  kinematic  calculations,  formation 
control  law  and  formation  separation  commands  are  implemented  in  the  wing  aircraft 
block/modtile.  Thus  additional  aircraft  can  be  added  to  the  formation  simulation  by 
duplicating  the  Wing  Aircraft  block.  This  is  explained  in  detail  in  Chapter  IV. 
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Figure  2.1.  Block  Diagram  of  Formation  Control  System  Simulation 
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The  following  steps  describe  in  detail  the  simulation  of  a  typical  formation  maneuver 
initiated  by  the  lead  aircraft: 

1.  The  pilot  of  the  lead  aircraft,  block  1,  commands  a  velocity,  heading,  or  altitude 
change.  The  lead  aircraft,  block  2,  responds  to  this  command. 

2.  The  sensors  on  the  wing  aircraft,  block  3,  relay  the  lead  aircraft’s  positional  data 
to  the  kinematics  calculations,  block  4.  This  block  calculates  the  actual  separation 
distances  and  the  lead  aircraft’s  position,  velocity,  and  heading.  This  information  is 
relayed  to  the  formation  controller. 

3.  The  formation  controller,  block  5,  receives  nominal  separation  information  from  block 
6,  and  kinematic  information  from  block  4.  An  error  signal  is  constructed  by  extract¬ 
ing  the  difference  between  the  commanded  and  actual  separation  parameters. 

4.  The  formation  controller  implements  a  Proportional  Plus  Integral  control  law,  com¬ 
manding  the  wing  aircraft,  block  7,  in  order  to  null  the  error  signal. 

5.  This  process  continues  as  long  as  the  formation  control  system  is  eng^ed. 

If  a  formation  separation  change  is  required,  then  block  6,  the  nominal  separation 
commands,  would  send  this  information  to  the  formation  controller.  The  following  sections 
describe  in  detail,  the  key  blocks  comprising  the  formation  control  system  simulation. 


2.1  Aircraft  and  Autopilot  Models 

2.1.1  First  Order  Models  A  major  component  of  the  formation  control  system 
simulation  are  the  aircraft  and  autopilot  models,  blocks  2  and  7  in  Figure  2.1.  Previous 
research  by  Dargan(2)  and  Ilohs(9)  developed  first-order  transfer  functions  modelling  the 
response  of  autopilot  equipped  C-130H  adrcraft.  The  lateral  and  longitudinal  autopilots 
were  divided  into  three  separate  channels,  velocity,  (magnetic)  heading,  and  altitude.  The 
transfer  functions  for  these  channels  are  shown  in  Eqs  (2.1)  to  (2.3)  (2)[3-3]. 


M.  = 


(2.1) 


2-2 


where 


V>(s)  _ 

_1_ 

M^) 

Ks)  _ 

Tk 

he{s) 

(2.2) 

(2.3) 


V  =  aircraft’s  velocity 

Ve  =  aircraft’s  velocity  command 

Tv  =  velocity  time  constant  of  first-order  aircraft  model 
0  =  aircraft’s  heading 

=  aircraft’s  heading  command 

=  heading  time  constant  of  first-order  aircraft  model 
h  =  aircraft’s  altitude 
he  =  aircraft’s  altitude  command 
T/,  =  altitude  time  constant  of  first-order  aircraft  model 

The  performance  variation  between  the  three  channels  dynamics  is  due  to  the  differences 
in  their  first-order  time  constants  tv,t^  and  respectively.  Additionally,  two  models 
of  the  C-130,  the  C-130A  and  the  C-130B,  are  simulated.  While  the  same  linear  transfer 
functions  are  used  for  both  aircraft,  they  differ  in  their  velocity,  heading  and  climb  rate 
limits(2)[3-2].  Table  2.1  displays  the  time  constants  and  rate  limits  for  each  aircraft. 

2.1.2  Second  Order  Models  Second  Order  Aircraft  models  are  based  on  second 
order  linear  transfer  functions  that  incorporate  a  time  delay  and  a  maximum  rate  of  change. 
These  transfer  functions  can  better  model  physical  systems  that  exhibit  an  appreciable 
amount  of  delay  in  their  time  response.  Additionally,  a  second-order  transfer  function  can 
model  both  an  underdamped  or  an  overdamped  time  response. 

A  major  difficulty  in  deriving  second  order  autopilot  models  of  the  C-130  aircraft 
is  the  acquisition  of  accurate  flight  performance  data.  Consequently,  the  second  order 
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Table  2.1.  Aircraft  Rate  Limits  and  Model  Time  Constants 


Aircraft 

and 

Lower 

Upper 

Tv 

Th 

Parameter 

Limit 

Limit 

C-130A 


Constants 


Velocity 


Heading 


Altitude 


C-130B 


Constants 


Velocity 


Heading 


Altitude 


.333  secs  1.5  secs 


304  ft/s 


-3  deg/s 


-42  ft/s 


422  ft/s 


3  deg/s 


8.5  ft/s 


.333  secs  1.5  secs  2  secs 


304  ft/s  422  ft/s 


-4.7  deg/s  4.7  deg/s 


-42  ft/s  33  ft/s 


transfer  functions  are  developed  from  computer  simulations  of  linear  bare  aircraft  state 
space  models,  equipped  with  basic  autopilots.  The  theoretical  method  for  synthesizing 
the  augmented  transfer  functions  is  contained  in  Appendix  A.  Appendix  B  describes  in 
detail,  the  autopilot  design  and  second  order  transfer  function  derivation.  The  new  transfer 
functions  for  the  C-130  axe  shown  below: 


II 
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where 

=  second  order  real  pole,  heading  channel 
=  second  order  real  pole,  heading  channel 
Oh  =  second  order  real  pole,  altitude  channel 
bh  =  second  order  real  pole,  altitude  channel 
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Of  special  note  is  the  velocity  chajinel  transfer  function,  Eq  (2.4).  This  channel  is  best 
modeled  by  a  first  order  transfer  function  similar  to  Eq  (2.1),  except  that  now  the  (real¬ 
istic)  time  constant  is  10  seconds,  as  opposed  to  0.333  seconds.  The  same  rate  limiters 
distinguishing  the  C-130A  from  the  C-130B  models  are  used  with  the  second  order  transfer 
functions.  A  comparison  of  the  new  and  previous  autopilot  models  time  responses  to  a 
step  input  is  shown  in  Figure  2.2.  The  time  responses  of  the  new  autopilot  models  are 


Figure  2.2.  C-130B  Aircraft  Autopilot  Models  Time  Response  Comparison 

slower  than  the  previous  models.  Consequently,  the  controller  parameters  determined  in 
previous  research  need  to  be  reevaluated.  The  first  order  models  are  used,  due  to  their 
simplicity,  for  theoretical  analysis  and  performance  comparison  against  previous  research. 
The  new  models  are  used  to  validate  the  final  control  system  parameters. 


2.2  Aircraft  Sensor  Measurements 

An  on-board  sensor  capable  of  providing  relative  position  and  velocity  information  of 
the  lead  aircraft  with  respect  to  the  wing  is  needed  (2)[3-8].  A  Kalman  Filter  can  be  used 
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to  process  the  available  measurements  in  order  to  estimate  non-measured  parameters  and 
improve  the  accuracy  of  the  measured  data.  For  the  purposes  of  the  computer  simulation, 
perfect  measurement  data  is  assumed. 

2.3  Formation  Coordinate  System 

Two  coordinate  reference  frames,  an  inertial  and  a  rotating  reference  frame  centered 
on  the  wing  aircraft,  are  used  in  this  research(2)[3-8]. 

The  inertial  reference  frame  has  a  stationary  origin,  with  latitude,  longitude,  and 
altitude  as  its  axes(2)[3-9].  During  computer  simulation,  the  lead  aircraft  begins  all  ma¬ 
neuvers  on  the  inertial  axis,  at  an  altitude  of  500  feet,  and  a  zero  degrees  magnetic  heading. 
This  somewhat  simplifies  the  kinematic  calculations.  Figure  2.3  depicts  the  lead  and  wing 
aircraft  separation  distances  using  the  inertial  frame  of  reference. 

The  rotating  reference  frame  origin  is  aligned  with  the  wing  aircraft,  with  the  x-axis 
aligned  to  the  aircraft’s  direction  of  flight,  y-axis  out  the  right  wing,  and  the  z-axis  down 
toward  the  earth(2)[3-9].  The  reference  frame  rotates  with  respect  to  the  inertial  reference 
frame,  as  the  wing  aircraft  changes  heading.  An  example  of  the  wing’s  rotating  reference 
frame  is  shown  in  Figure  2.4.  The  orthogonal  separation  distances  between  the  wing  and 
lead  aircraft  are  directly  provided  by  this  coordinate  system. 

Both  reference  frames  are  used  extensively  throughout  this  research.  Kinematic 
equations  are  developed  in  the  next  section  utilizing  these  coordinate  systems. 

2.4  Kinematic  Equations 

The  kinematic  equations  are  used  to  represent  the  relative  positions/distances  be¬ 
tween  the  lead  and  wing  aircraft.  The  distances  are  expressed  using  a  Cartesian  coordinate 
system  affixed  to  the  inertial  reference  frame  or  the  wing’s  rotating  reference  frame.  In  pre¬ 
vious  research  by  Dargan,  the  separation  distances  between  the  two  aircraft  with  respect 
to  the  wing  was  derived(2)[3-12  to  3-21].  The  lead-wing  separation  distances,  referred  to  as 
the  inertial  reference  frame,  is  derived  in  Chapter  IV.  The  lead-wing  separation  distances, 
referred  to  as  the  wing  aircraft,  is  reviewed  in  the  following  paragraphs.  Given  the  inertial 
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Y  Inertial 


Figure  2.3.  Inertial  Reference  Frame  and  Separation  Distances 


velocity,  heading  and  altitude  of  each  aircraft,  an  equation  which  relates  the  position  of 
the  lead  aircraft  with  respect  to  the  wing  is  desired.  Using  the  Equation  of  Coriolis,  the 
velocity  of  lead  aircraft  with  respect  to  the  wing  is  shown  below(2) [Equation  (3.12)]: 


^WL 


X  Rwl 


_  V"' 
yvL 


(2.7) 


where 


V^i  =  lead’s  velocity  with  respect  to  the  wing,  expressed  in  wing’s  frame 
=  lead’s  inertial  velocity  in  wing’s  frame 
=  angular  velocity  of  wing’s  frame  expressed  in  wing’s  frame 
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Figure  2.4.  Wing’s  Rotating  Reference  Frame  and  Separation  Distances 


Rwi  =  position  of  lead  with  respect  to  wing  expressed  in  wing’s  frame 
—  wing’s  velocity  expressed  in  its  own  frame  of  reference 


and, 


where 


VH'  =  c^vt  = 


cos  V’fi  —  sin  ipE  0 
sin  xpE  cos  ipE  0 
0  0  1 


Ve 

Vl  cos  XpE 

0 

= 

Vl  sin  xpE 

0 

0 

=  direction  cosine  matrix,  lead  to  wing  transformation 
tpE  =  ‘•Pl-  i’w,  heading  error  between  lead  and  wing 


(2.8) 
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Y  Inertial 


Separating  the  above  equation  into  the  y,  and  z  components  yields: 


=  Vl  cos IpE  +  ^wy'^  -  Vw 

(2.11) 

=  Vi  sin  ipB  - 

(2.12) 

=  0 

(2.13) 

In  the  computer  simulation,  the  above  equations  are  integrated  to  get  the  separation 
distances  as  a  function  of  time. 

2.5  Formation  Controller  Strategy 

In  this  research,  a  two  tiered  control  strategy  is  employed(2)[3-ll].  The  upper  tier 
is  the  control  of  the  formation  as  a  whole,  and  the  lower  tier  is  for  the  control  of  the 
individual  aircraft  within  the  formation.  In  Figure  2.1,  the  upper  tier  formation  commands 
are  directed  by  the  pilot  in  the  lead  aircraft,  while  the  lower  tier  commands  are  directed  by 
the  formation  controller  located  on  each  wing  aircraft.  This  section  describes  the  design 
of  the  formation  controller. 

In  previous  research  by  Dargan,  two  types  of  formation  controllers  utilizing  propor¬ 
tional  plus  integral  feedback  were  designed.  The  first  type  utilized  only  a  PI  controller, 
while  the  second  type  included  both  a  PI  controller  and  a  linear  mixer.  Both  types  of 
formation  controllers  are  used  in  this  research. 

2.5.1  Closed  Loop  Control  Using  a  PI  Controller  The  PI  controller  operates 
on  the  error  between  the  commanded  separation  and  the  actual  separation  distances.  The 
proportional  and  integral  action  is  used  to  drive  this  error  to  zero.  A  block  diagram  of  the 
Formation  Control  System  utilizing  the  PI  Controller  is  shown  in  Figure  2.6. 
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where 


rl  = 


Lead  Velocity 
Lead  Heading 
Lead  Altitude 


yl  = 


Wing  Velocity 
Wing  Heading 
Wing  Altitude 


r2  = 


X  Separation 
Y  Separation 
Z  Separation 


2/2  = 


Wing  X  Position 
Wing  Y  Position 
Wing  Z  Position 


(2.14) 


2/3  = 


Wing  Velocity  Command 
Wing  Heading  Command 
Wing  Altitude  Command 


el  =  rl  —  yl  =  Maneuver  Errors 
e2  =  r2  —  y2  =  Separation  Errors 


A  large  command  is  supplied  to  the  wing  aircraft  when  either  the  magnitude  or  inte¬ 
gral  of  the  error  is  large(2)[4-l].  The  key  to  this  controller’s  operation  is  the  proportional 
plus  integral  gains  in  the  PI  controller.  The  following  equations  define  the  relationship 
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between  the  lon^tudinal  and  lateral  separation  errors  and  the  command  inputs  to  the 
wing  aircraft. 


Pive(0  —  kgpX  -}-  k^i  1  xdt 

Jo 

(2.15) 

ipwcit)  =  kypy  -1-  kyi  /  ydt 

Jo 

(2.16) 

hn^^(t)  =  kipZ  -f-  kii  I  zdt 

Jo 

(2.17) 

(2.18) 

where 

X  —  —  x^ ,  X  separation  error 

y  =  —  y^  ,y  separation  error 

2  =  ,z  separation  error 

kgp  =  x-channel  proportional  gain 
kxi  =  x-channel  integral  gain 
kyp  =  y-channel  proportional  gain 
kyi  =  y-channel  integral  gain 
kip  =  z-channel  proportional  gain 
kii  =  z-channel  integral  gain 

The  proportional  and  integral  controller  gains  are  optimized  in  order  to  produce  the  best 
time  response.  The  previous  research  by  Dargan  showed  that  this  type  of  formation  control 
system  produced  adequate  transient  responses  and  zero  steady  state  error  (2)[4-13]. 

2.5.2  Closed  Loop  Control  Using  a  Linear  Mixer  and  a  PI  Controller  The 
following  formation  control  system  is  similar  to  the  previous  one,  except  a  linear  mixer  is 
added.  The  mixer  is  used  to  combine  the  separation  errors  (x  and  y  spacing  errors)  as  well 
as  the  maneuver  errors  (velocity  and  heading  errors).  The  PI  controller  operates  on  this 
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combined  error,  by  commanding  the  wing  aircraft  to  drive  this  error  to  zero.  The  addition 
of  the  velocity  and  heading  error  signal  improves  the  transient  response  of  the  system.  The 
key  parameters  to  the  controller  operation  are  the  gains  used  in  both  the  mixer  and  the  PI 
controller.  Only  velocity  and  heading  channels  are  controlled  using  this  formation  control 
system.  A  block  diagram  of  the  Formation  Control  System  utilizing  the  mixer  and  the  PI 
Controller  is  shown  in  Figure  2.7. 


where 


rl  = 


Lead  Velocity 
Lead  Heading 


r2  = 


X  Separation 
Y  Separation 


y\  = 


Wing  Velocity 
Wing  Heading 


y2  = 


Wing  X  Position 
Wing  Y  Position 


(2.19) 


yZ  = 


Wing  Velocity  Command 
Wing  Heading  Command 


el  =  rl  —  yl  =  Maneuver  Errors 
e2  =  r2  -  y2  =  Separation  Errors 
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The  following  equations  dehne  the  relationship  between  the  lon^tudinal  and  lateral 
separation  errors  and  the  command  inputs  to  the  wing  aircraft. 

^We(0  ~  +  kxx]  +  k,;i  j  [AivVe  +  kxx]  dt  (2.20) 

Jo 

ipwcit)  =  kyp[k^xl)E  +  kyy]  +  kyi  f  +  fcyy]  dt  (2.21) 

Jo 

where 

Ve  =  Vi,-  Fw,  velocity  error  signal 
kv  =  velocity  error  signal  mixer  gain 
k^  =  heading  error  signal  mixer  gain 
kx  =  x-separation  error  signal  mixer  gain 
ky  =  y-separation  error  signal  mixer  gain 

The  mixer  and  PI  controller  gains  are  again  optimized  in  order  to  produce  the  best  time 
response.  The  use  of  the  mixer  provides  a  better  transient  response  behavior  for  the 
longitudinal  channel  but  there  is  little  improvement  in  the  lateral  channel  (2)[4-22]. 

2. 6  Measures  of  Merit 

The  measures  of  merit  are  based  on  the  formation  control  system’s  transient  behavior 
and  steady  state  error  performance.  The  formation  control  system  is  considered  to  perform 
satisfactorily  if  it  is  able  to  track  all  commanded  inputs  with  zero  steady  state  error,  and 
if  its  transient  behavior  is  such  that  “collisions”  between  the  lead  and  the  wing  aircraft 
are  avoided  during  the  transients(2)[3-36].  An  important  consideration  is  the  stability  of 
the  control  system.  Large  formation  maneuver  command  inputs,  such  as  a  90  degree  step 
heading  change,  can  produce  unstable  responses  for  large  value”  of  controller  gains,  due  to 
the  deleterious  effects  of  the  system’s  nonlinearities.  The  ability  of  the  formation  control 
system  to  handle  large  command  inputs  is  also  used  as  a  figure  of  merit. 
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III.  Lead  Aircraft’s  TYajectory  in  Wing  Aircraft’s  Rotating  Reference 

Frame 


The  performauice  of  the  formation  control  system  is  determined  by  analyzing  the 
time  responses  of  the  wing  aircraft.  The  wing  aircraft’s  velocity  and  heading  formation 
separation  distances,  as  well  as  altitude  is  monitored.  These  distances  are  of  special  concern 
since  they  determine  how  well  the  formation  is  maintained  during  maneuvers  and  how 
well  formation  geometry  change  commands  are  executed.  While  the  time  response  plots 
provide  valuable  information  on  the  transient  and  steady  state  behavior  of  the  formation 
control  system,  they  fail  to  provide  insight  into  the  spatial  trajectory  of  the  wing  aircraft. 
Presenting  the  wing  or  lead  aircraft’s  trajectory  with  respect  to  the  other  aircraft  provides 
a  visualization  tool  in  evaluating  the  performance  of  the  formation  control  system. 

3.1  Lead  Aircraft’s  Trajectory  using  Lissajous  Figures 

The  X  and  y  separation  distances  between  the  lead  and  wing  aircraft  are  obtained  as 
a  function  of  time,  with  respect  to  the  wing’s  rotating  reference  frame.  These  distances  are 
presented  in  Eqs  (2.11)  and  (2.12).  Using  this  information,  the  lead  aircraft’s  trajectory, 
which  is  given  in  parametric  form,  is  plotted  using  a  Cartesian  coordinate  system.  Figure 

3.1  illustrates  a  two  dimensional  trajectory.  The  wing  aircraft  is  located  at  the  origin  of 
the  diagram.  The  solid  line  in  Figure  3.1  is  the  trajectory  of  the  lead  aircraft  in  the  wing 
aircraft’s  reference  frame.  The  trajectory  is  a  closed  path  since  the  original  formation  is 
maintained  after  the  maneuver  is  complete.  This  would  not  be  the  case  during  a  formation 
change  maneuver.  The  dotted  circle  represents  a  minimum  formation  separation  distance 
which  is  established  as  a  safety  region.  If  the  solid  line  were  to  penetrate  this  circle,  then 
the  minimum  separation  distance  specification  is  violated. 

This  presentation  of  the  lead  aircraft’s  trajectory  is  analogous  to  Lissajous  figures 
used  to  visualize  two-dimensional  oscillations.  The  shape  of  this  figure  provides  insight 
in  the  response  of  the  formation  control  system.  A  small  figure  indicates  a  tight  and  fast 
control  system  with  little  deviation  from  the  original  formation.  A  large  and  wide  figure 
indicates  a  sluggish  response  with  wide  swings  in  the  formation. 
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y  separation  (feet) 

Figure  3.1.  Example  of  Lead  Aircraft’s  Trajectory  Lissajous  Figure 

3.2  Performance  Evaluation  using  Lissajous  Figures 

3.2.1  Diamond  Formation  Heading  Change  The  time  response  plots  of  a  sin¬ 
gle  wing  aircraft  formation  through  a  plus  and  minus  30  degree  heading  change  maneuver 
is  shown  in  Fiqure  3.2.  Two  separate  simulations  are  included  in  the  figure.  The  solid 
lines  represent  the  plus  30  degree  heading  change,  while  the  dotted  lines  are  the  negative 
30  degree  heading  change.  This  information  suffices  to  simulate  a  symmetric  three  aircraft 
diamond  formation’s  heading  change.  The  corresponding  Lissajous  figures  are  shown  in 
Figure  3.3.  In  both  maneuvers,  the  lead  mrcraft’s  trajectory  is  a  closed  path  since  the  for¬ 
mation  is  maintained  after  the  maneuver.  The  dotted  line  represents  a  600  feet  separation 
boundary  between  the  two  aircraft.  The  lead  aircraft  is  commanded  inward  toward  the 


3-2 


wing  during  the  minus  30  degree  heading  maneuver,  consequently  its  trajectory  approaches 
closer  to  the  600  ft  boundary. 


Figure  3.2.  Diamond  Formation,  ±30  Degrees  Heading  Change  Time  Response 
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separation  (feet)  x  seporotion  (feet) 
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y  seporotion  (foot) 

Figure  3.3.  Diamond  Formation,  +30  Degrees  Heading  Change  (top),  and  —30  Degrees 
Heading  Change  (bottom).  Lead  Trajectory  Lissajous  Figures 
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3.2.2  Trail  Formation  Heading  Change  The  previous  simulations  are  repeated 
for  two  aircraft  in  a  trail  formation.  The  plus  and  minus  30  degree  heading  change  time 
response  plots  of  the  two  aircraft  are  shown  in  Figure  3.4.  The  lead  aircraft’s  trajectories 
are  shown  in  Figure  3.5.  As  for  the  case  for  the  diamond  formation  heading  change, 
the  lead  aircraft’s  trajectory  is  a  closed  path,  since  the  formation  is  maintained  after  the 
maneuver.  A  500  feet  boundary  is  established  for  this  maneuver.  As  seen  in  Figure  3.5,  the 
boundary  is  crossed  during  both  heading  changes.  Even  though  the  boundary  is  violated, 
the  magnitude  of  the  crossing  is  less  than  one  foot.  As  expected,  the  Lissajous  figures  for 
both  are  identical,  but  are  reflected  about  the  x  axis. 


0  10  20  30  40  so  60 

Figure  3.4.  Trail  Formation, ±30  Degrees  Heading  Change  Time  Response 
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3.3  Summary 

The  Lissajous  figure  representation  of  the  lead  aircraft’s  trajectory  in  the  wing  air¬ 
craft’s  rotating  reference  frame  is  a  valuable  tool  for  visualizing  the  performance  of  the 
formation  control  system.  It  provides  a  graphical  means  of  verifying  any  violation  of  the 
minimum  distance  between  the  two  aircraft.  The  size  and  shape  of  the  lead  aircraft’s  tra¬ 
jectory  is  used  to  evaluate  the  tightness  of  the  formation  and  speed  of  the  control  system. 
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IV.  Two  Wing  Aircraft  Formation  Flight  Control 


A  degree  of  complexity  is  added  to  the  formation  flight  control  problem  for  formations 
including  two  wing  aircraft.  The  two  wing  aircraft  not  only  must  independently  and 
successfully  track  the  lead  aircraft,  but  they  must  also  avoid  a  collision  between  themselves. 
The  current  formation  control  concept  stipulates  that  each  wing  aircraft,  through  the  use 
of  its  sensors,  tracks  the  lead  aircraft  by  measuring  the  lead-wing  separation  distances. 
The  separation  distance  between  the  two  wing  aircraft  is  not  measured,  but  it  is  computed 
analytically  for  this  report.  Furthermore,  no  communication  nor  coordination  between  the 
wing  aircraft  is  envisioned.  However,  the  formation  control  simulation  is  used  to  ensure  no 
collision  occurs  between  the  two  wing  aircraft  during  a  formation  maneuver  or  a  formation 
geometry  change.  This  is  the  primary  motivation  for  the  implementation  of  the  two  wing 
formation  flight  control  system  simulation. 

4.1  Two  Wing  Formation  Controller  Implementation 

The  Two  Wing  Formation  Flight  Control  Problem  is  simplified  through  the  use  of 
the  following  assumptions. 

•  Each  wing  aircraft  is  equipped  with  a  formation  flight  control  system  capable  of 
independently  tracking  the  lead  aircraft  with  an  established  tracking  error. 

•  The  separation  distances  determining  the  geometry  of  the  formation  are  established 
in  accordance  with  the  tracking  error  of  each  wing  aircraft.  For  example,  a  heavier 
aircraft  with  a  high  tracking  error  requires  larger  separation  distances. 

•  These  assumptions  can  be  used  in  establishing  multiple  wing  aircraft  formations. 

The  diagram  in  Figure  4.1  outlines  the  geometry  of  a  two  wing  aircraft  formation.  In 
this  diagram,  a  diamond  formation  is  illustrated.  The  addition  of  a  second  wing  aircraft 
requires  the  use  of  a  common  inertial  reference  frame,  and  the  calculation  of  the  inertial 
positions  of  each  aircraft.  This  common  reference  frame  is  used  to  obtain  the  separation 
distances  between  the  three  aircraft  during  simulations.  The  calculations  are  presented  in 
the  next  section. 
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y  Inertial 

Figure  4.1.  Two  Wing  Aircraft  Formation  Geometry 


4.S  Inertial  Calculations 

The  inertial  position  of  each  aircraft  is  needed  in  order  to  calculate  the  separation 
distance  between  two  aircraft.  Depending  on  the  information  available,  the  inertial  sepa¬ 
ration  distance  between  aircraft  can  be  calculated  in  two  different  ways.  In  the  computer 
simulation,  the  inertial  velocity  of  each  aircraft,  resolved  in  its  own  reference  frame,  and 
the  inertial  headings  are  known.  AdditionaJly,  the  initial  inertial  position  of  each  aircraft 
is  known.  Without  loss  of  generality,  the  initial  position  of  the  lead  aircraft  is  assumed 
to  be  the  inertial  origin.  Equations  (4.2)  to  (4.6)  generate  the  positions  of  each  aircraft. 
The  separation  distances  between  the  aircraft  is  the  difference  between  these  values.  Since 
the  axis  of  each  aircraft  is  aligned  with  the  axis,  no  conversion  is  necessary.  The 
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inertial  separation  between  the  two  wing  aircraft  is  shown  in  Eq  (4.9). 


XL  =  f  Vicos i>i,  +  Xi(0) 

Jto 

Yl  =  /V/'sin^i.  +  r/CO) 

Jto 

Xwi  —  f  V^' cosipwi  + 

Jto 

Yw.  =  /V“;‘sin^v„.  +  y^,(0) 

Jto 

Xw,  =  I  V^^  cos  X^fJ^Qi) 

Jto 

Ywo  =  sin  + 5-4(0) 

Jto 

where 


XL  =  lead’s  inertial  position  along  X^  axis 
¥[  =  lead’s  inertial  position  along  axis 
Xj/Vi  —  wing'jS  inertial  position  along  axis 
5^w,  =  wingjS  inertial  position  along  Y^  axis 
V’w,  =  wingJs  inertial  heading 
Xwo  =  wingjS  inertial  position  along  Y^  axis 
5^4  =  wingjS  inertial  position  along  Y^  axis 
i’Wi  =  wingjS  inertial  heading 
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(4.1) 

(4.2) 

(4.3) 

(4.4) 

(4.5) 

(4.6) 


(4.7) 


(4.8) 

(4.9) 


where 


=  lead-wingi  inertial  separation 
=  lead-wing2  inertial  separation 
RwiWa  =  wing2-wingi  inertial  separation 


An  alternative  method  of  calculating  the  inertial  separations  between  the  aircraft  uses 
the  separation  distance  between  the  wing  and  the  lead  aircraft  referenced  to  the  wing’s 
rotating  reference  frame.  In  an  actual  implementation,  this  information  is  available  to  the 
each  wing  aircraft.  The  inertial  separation  between  the  two  aircraft  is  calculated  using 
the  direction  cosine  matrix  in  Eq  (4.10).  Equations  (4.11)  to  (4.12)  implement  the  steps 
needed  to  calculate  the  inertial  separations. 


cos  V’w. 

-  sin  V’w. 

0 

a:n 

— 

- 

sinV»w. 

cos 

0 

R*^-  - 

l»n 

0 

0 

1 

(4.10) 


where 


=  direction  cosine  matrix,  nth  wing  to  inertial  frame  of  reference  (f.o.r.) 
=  nth  wing  aircraft’s  inertial  heading 
x„  =  nth  wing  x-separation  in  its  f.o.r. 
yn  =  nth  wing  y-separation  in  its  f.o.r. 

Zn  =  nth  wing  z-separation  in  its  f.o.r. 

R^“i  =  lead’s  position  with  respect  the  nth  wing  aircraft 


and 


iicos^vv,  -  l/i  sinV’w, 
xi  sin  V’w,  +  yi  cos  V’w, 


(4.11) 
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X2  COS  rpw,  -  Vi  sin  ipw, 

—  X2Sinxj>W3  +  I/l  cos  V’lVa  (4-12) 

22 

The  inertial  separation  between  the  two  wing  aircraft  is  the  same  as  Eq  (4.9).  The 
distance  can  be  transformed  to  wing(s  reference  frame  using  the  direction  cosine  matrix. 
The  transformation  from  the  inertial  reference  frame  to  wingjs  reference  frame  is  shown 
in  Eq  (4.13). 

(4-13) 

This  distance  is  useful  in  plotting  wingjS  trajectory  with  respect  to  wingi.  A  Lissajous 
figure  showing  this  trajectory  is  included  in  the  next  section. 

4.3  Performance  Evaluations 

The  lead  aircraft  is  commanded  a  30  degree  heading  change  while  the  wing  aircraft 
are  in  both  trail  and  diamond  formations.  During  these  simulations,  no  altitude  change 
is  commanded,  so  the  z  coordinate  is  ignored.  The  separation  distance  between  all  three 
aircraft  is  observed  in  order  to  meet  minimum  separation  requirements. 

4-3.1  Trail  Formation  Heading  Change  The  time  response  plots  of  a  two  wing 
aircraft  trail  formation  through  a  30  degree  heading  change  maneuver  is  shown  in  Fig¬ 
ure  4.2.  All  aircraft  are  initially  separated  by  500  feet.  Wingi  is  directly  behind  the  lead 
aircraft,  and  wingj  is  500  feet  behind  wingi.  See  Fiqure  1.1  for  reference.  The  time  re¬ 
sponses  of  both  wing  aircraft  are  included  in  this  figure.  The  separation  distances  shown 
in  Figure  4.2  are  with  respect  to  each  wing  aircraft. 

The  inertially  resolved  separation  distances  between  all  three  aircraft  during  the  30 
degree  heading  change  are  shown  in  Figure  4.3.  The  closest  distance  between  any  two 
aircraft,  for  this  maneuver,  is  between  the  wingi  and  wingo.  The  separation  distance 
between  these  two  aircraft  decreases  from  a  nominal  500  ft  to  494  ft.  Even  though  there 
is  no  direct  control  of  this  parameter,  the  individual  formation  control  system  on  each 
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tinM(aac) 


Figure  4,2.  Trail  Formation,  Two  Wing  Aircraft,  +30  Degrees  Heading  Change  Time 
Response 


aircraft  is  able  to  maintain  the  commanded  separation.  Figure  4.4  displays  a  Lissajous 
figure  containing  wingjS  trajectory  in  wingjS  reference  frame. 
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Figure  4.3.  Trail  Formation,  Two  Wing  Aircraft,  +30  Degrees  Heading  Change,  Separa¬ 
tion  Distances 


Figure  4.4.  Lissajous  Figure  of  WingjS  lYajectory  in  Wingjs  Reference  Frame,  for  a  +30 
Degrees  Heading  Change  in  Trail  Formation 
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4.3.2  Diamond  Formation  Heading  Change  The  tin  e  i?sponse  plots  of  a  two 
wing  aircraft  diamond  formation  through  a  30  degree  heading  cha^ige  maneuver  is  shown 
in  Figure  4.5.  The  radial  separation  distances  is  shown  in  Figure  4.6.  Of  special  note  is 
the  separation  of  the  two  wing  aircraft.  The  minimum  distance  between  the  two  aircraft 
is  just  under  996  ft.  There  is  little  deviation  from  the  nominal  1000  feet  separation  during 
this  maneuver. 


-10  O  10  20  30  40  50  60  70  80 

y1  -2  (feet) 

Figure  4.5.  Diamond  Formation,  Two  Wing  Aircraft,  +30  Degrees  Heading  Change  Time 
Response 

The  relatively  high  tracking  performance  of  the  formation  controllers  on  each  aircraft 
is  able  to  maintain  the  distance  between  all  aircraft.  This  is  also  illustrated  in  Figure  4.7. 
In  this  figure,  a  Lissajous  figure  containing  wingjS  trajectory  with  respect  to  wingi  is 
displayed.  The  next  plot.  Figure  4.8,  displays  the  flight  path  of  each  aircraft  during  the  30 
degree  heading  change.  This  inertial  perspective  illustrates  how  each  wing  aircraft  is  able 
to  track  the  lead.  The  small  circles  represent  five  second  intervals.  They  can  be  used  to 
roughly  show  the  changing  position  of  each  aircraft. 
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Figure  4.6.  Trail  Formation,  Two  Wing  Aircraft,  +30  Degrees  Heading  Change,  Separa¬ 
tion  Distances 


975  980  985  990  995  1000  1005 

y1  -2  (feet) 

Figure  4.7.  Lissajous  Figure  of  WingjS  Trajectory  in  Wingjs  Reference  Frame,  for  a  +30 
Degree  Heading  Change  in  Diamond  Formation 
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Figure  4.8.  Diamond  Formation,  +30  Degrees  Heading  Change,  Inertial  Flight  Paths 
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4-4  Summary 

In  this  chapter,  a  second  wing  aircraft  is  added  to  the  flight  control  problem.  This 
second  aircraft  is  equipped  with  a  formation  flight  controller  identical  to  the  first  wing 
aircraft.  In  order  to  properly  maintain  the  formation  of  the  two  wing  aircraft,  the  separa¬ 
tion  between  these  aircraft  must  be  maintained.  Even  though  this  separation  distance  is 
not  directly  measured  and  controlled,  the  ability  of  each  wing  aircraft  to  successfully  track 
the  leader,  is  sufficient  to  maintain  the  formation,  and  guarantee  the  separation  of  wingi 
and  wing2.  The  wing  tracking  errors  due  to  the  heading  change  are  small  in  magnitude 
compared  to  the  formation  geometry.  This  is  verified  by  the  30  degree  heading  change 
computer  simulations.  Consequently,  the  formation  flight  control  system  developed  for  the 
single  wing  aircraft  can  successfully  be  used  for  midti-wing  aircraft  formations.  This  is 
true  as  long  as  the  individual  tracking  errors  of  each  aircraft  are  small  compared  to  the 
size  of  the  formation. 
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V.  Analytical  Formation  Flight  Controller  Synthesis 


The  formation  flight  controller  entails  output  feedback  and  proportional  plus  integral 
control.  This  combination  produces  output  decoupling  and  tracking  of  the  command  inputs 
by  the  corresponding  outputs.  Unfortunately  the  non-linear  rate  limiters  in  the  aircraft 
models  and  trigometric  functions  in  the  kinematic  equations  preclude  the  exclusive  use 
of  linear  control  design  methods.  Hence,  both  the  proportional  and  integral  gains  in  the 
PI  controller  are  determined  using  the  formation  flight  control  simulation  and  a  trial  and 
error  method.  Numerous  computer  simulations  are  performed  while  modifying  these  gains 
and  observing  the  tracking  ability  of  the  control  system.  This  method  produces  successful 
results,  but  considerable  labor  is  expended  in  determining  the  PI  controller  gains. 

A  linearized  plant  model  is  desirable  since  pole  placement  techniques  can  be  used  to 
analytically  determine  initial  controller  gain  values.  These  initial  values  are  then  experi¬ 
mentally  refined  using  the  full  non-linear  computer  simulation.  This  chapter  presents  the 
derivation  of  a  linear  and  parametrized  MIMO  plant  used  to  analytically  determine  the 
optimal  PI  control  gains  for  the  formation  flight  control  problem.  In  this  chapter,  analysis 
is  restricted  to  the  planar  velocity  and  heading  channels. 

5.1  Theoretical  Development 

The  linear  model  for  the  formation  flight  control  problem  is  developed  as  follows: 
First,  the  formation  geometry  is  outlined  and  the  key  geometric  parameters  are  identified. 
Linearized  equations  for  both  the  aircraft  models  and  kinematics  are  developed.  These 
equations  are  then  non-dimensionalized  and  presented  in  state  space  form.  The  separate  x 
and  y-channels  are  identified,  decoupled,  and  augmented  using  the  PI  control  law.  These 
augmented  state  equations  are  then  analyzed  for  steady  state  error  and  stability.  Finally, 
the  characteristic  equation  of  these  augmented  plants  are  used  in  determining  the  optimal 
control  laws  (gains).  The  formation  geometry  is  illustrated  in  Figure  5.1.  The  key  param¬ 
eters  are  the  separation  I  =  y/x^  -f  y^,  and  the  formation  angle  a.  Equations  (5.1)  through 


5-1 


Figure  5.1.  Formation  Geometry 


(5.2)  relate  the  Cartesian  separation  distances  to  these  two  parameters, 


Xo  =  I  cos  a 
jfo  =  Isina 


where 


Xg  =  nominal  x-separation 
Vo  =  nominal  y-separation 
a  =  separation  angle 


First  order  aircraft  models  are  used  in  this  analysis.  This  simplifies  the  final  equations 
and  provides  a  means  of  comparing  the  final  results  with  the  resxilts  developed  by  Dargan 
(2)[4- 1-4-22]  and  Rohs  (9).  Equations  (5.3)  to  (5.5)  define  the  first  order  C-130  aircraft 
models  developed  in  Rohs  research  (9)[28-35]. 
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;  1  .  1  , 

"■nr  — - hw  H - nwe 

Th  W  Th  w 

where 

Tvyt,  =  wing’s  first  order  velocity  time  constant 
r^n,  =  wing’s  first  order  heading  time  constant 

The  kinematic  equations  describing  the  separation  between  the  lead  and  the  wing  aircraft 
are  given  by(2)[Eqs.  3.15  and  3.16] 


X  =  Vi  COSINE  +  i’wy  -  Vw 

(5.6) 

y  =  Vl  sin -ipE  -'tj’wx 

(5.7) 

Incorporating  small  angle  approximations  and  the  small  perturbations  method  into  Eqs.  (5.6) 

and  (5.7)  yields 

X  =  Vi,  +  ipwy„  -  Vw 

(5.8) 

y  =  Vl'^e  -  'I’wXo 

(5.9) 

Substituting  Eqs.  (5.1)  through  (5.4)  into  Eqs  (5.8)  and  (5.9)  yields 

X  =  V[,+ 

if  =  Vl  (V’j 

- —^w  +  — ^wc]  f  sin  Q  -  Vw 

[,  -  'Pw)  - - —i’W  +  -^V’We]  /  COS  a  -  vw 

(5.10) 

(5.11) 

(5.3) 

(5.4) 

(5.5) 
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The  next  step  in  the  development  of  the  linear  MIMO  plant  model  is  the  parametrization 
(non-dimensionalization)  of  all  variables.  Define  the  following  variables 


where 


(5.12) 

(5.13) 


/  =  characteristic  length,  nominal  formation  separation 

i  =  characteristic  time 

Vo  =  nominal  formation  velocity 

These  variables  are  used  to  non-dimensionalize  the  Eqs  (5.10)  and  (5.11).  Non-dimensional 
variables  and  parameters  are  defined  using  the  characteristic  length  and  time.  These 
variables  and  parameters  are  identified  using  the ‘symbol  and  are  displayed  below. 

^  -  I  y  -  I 

j)o  =  ^  =  sina  VW  =  ^ 

f  —  lYX.  f.  —  I*M. 

'Vw  —  f  >^w  ~  f 

where 

X  =  non-dimensional  x-separation  error 
y  =  non-dimensional  y-separation  error 
Xo  =  non-dimensional  nominal  x-separation 
yo  =  non-dimensional  nominal  y-separation 
Vyy  =  non-dimensional  wing  velocity 
Vwe  =  non-dimensional  wing  velocity  command 
fvyy  =  non-dimensional  wing  velocity  time  constant 
=  non-dimensional  wing  heading  time  constant 


Xo  =  ^  =  cos  a 

Vwc  =  ^  (5.14) 
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Vjr,  =  non-dimensional  lead  velocity 


The  above  non-dimensional  parameters  are  substituted  into  Equations  (5.3),  (5.4),  (5.10) 
and  (5.11)  n  order  to  form  the  following  linear  non-dimensional  differential  equations. 
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(5.15) 

(5.16) 

(5.17) 

(5.18) 


The  above  differential  equations  are  rearranged  into  state  variable  form.  Thus,  the  MIMO 
plant  is  represented  by  the  following  equations. 


X  =  AX  +  BU  +  TD 
Y  =  CX 


(5.19) 

(5.20) 


where  the  state  and  output  matrices  are 
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(5.21) 

1  0  0 

0 

Vc 

vw. 

X 

c  = 

,  £>  = 

,  u  = 

,  y  = 

(5.22) 

O 
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O 

_ 1 

0 

i’w 

C 

y 

where 


D  =  disturbance  input  matrix 
r  =  disturbance  input  vector 


5-5 


In  the  next  two  sections,  the  y  and  x-channels  are  identified  and  decoupled.  The  PI  control 
law  is  applied  to  these  decoupled  state  equations  to  form  the  closed  loop  state  equations. 

5.1.1  Y-Channel  The  y-channel  is  decoupled  from  the  x-channel  by  partitioning 
the  X  G  state  vector  as  shown  in  Eq  (5.23). 


X  =  [(i,VW),(j/,V»w)] 

Hence,  the  decoupled  y-channel  state  equation  is: 


(5.23) 


y 


0 

0 


COS  Qt  _ 

_ L_ 


*1 

■  * 

_ cos  Or 

1 

y 

+ 

^We  + 

1 

1 

0 

(5.24) 


The  wing  aircraft’s  heading  is  commanded  in  response  to  a  formation  separation  change 
in  the  y  direction.  The  PI  control  law  operates  on  the  lateral  separation  error  y  and  is 


■^wc  =  KpV  + 


/  ydt 
Jo 


(5.25) 


where 

kyp  =  non-dimensional,  y-separation  error,  proportional  gain 
kyi  =  non-dimensional,  y-separation  error,  integral  gain 


Differentiating  Eq  (5.25)  yields 

‘^Wc  —  kyiy  "h  ^ypy 
and  inserting  in  Eq  (5.26)  the  expression  for  y,  yields 


i’Wc  =  kyiV  + 


/cos  a 


*  cos  a 

V’W  -  kyp- - V’Wc  +  kyptpL 

T^yy 


The  augmented  state  equations  incorporating  the  PI  controller  are 


(5.26) 


(5.27) 


Xy  =  AyXy  +  Py  rf  j 


(5.28) 
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(5.29) 


i.e.. 


=  CyXy 


V’w 

= 

0 

0 


k  ■ 

n-yl 


coaor  _ 

^*w 

_ 1_ 

^*w 

w  (e  -  0 


cos  or 
**w 
1 


-k, 


yPr 


J  L 


y 

1 

+ 

0 

V’Wc 

1 

_ 1 

i’L 


(5.30) 


=  [» 


1  0 


y 

i’w 

i’We 


=  [i’w] 


(5.31) 


where 


Xy  =  augmented  y-channel  state  vector 
Ay  =  augmented  y-channel  state  transition  matrix 
Fj  =  augmented  y-channel  disturbance  input  matrix 
Dy  =  augmented  y-channel  disturbance  input  vector 
Yy  =  augmented  y-channel  output  vector 
Cy  =  augmented  y-channel  output  matrix 


Note:  The  two  y-channel  parameters  in  Eqs  (5.30)  and  (5.31)  are  the  formation  angle  a, 
and  the  heading  time  constant 

(i)  Statics  Consider  a  lead  heading  change  (disturbance)  input,  V’l(0-  The  steady 
state  response  to  a  unit  step  disturbance  is 


y{i).,  =  lims(5/-  ^)-‘r„i  =  -A7^ly 

5  9  • 


0 

1 

1 


(5.32) 
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The  wing  aircraft  must  be  able  to  track  the  lead  aircraft  with  zero  steady  state  error.  This 
zero  steady  state  tracking  error  requirement  is  stated  below, 


At  i  ->  oo,  y(oo)  =  0,  V’w{oo)  =  V’wcCoo)  =  =  Dy  =  1,  (5.33) 


i.e.. 


(5.34) 


Now  Eq  (5.33)  implies  that  Eq  (5.34)  satisfies  this  requirement,  as  long  as  Ay  is  invertible. 
A  necessary  condition  for  the  augmented  state  matrix  Ay  to  be  invertible  is  that  the  integral 
gain  kyi  must  be  non-zero. 


(ii)  Dynamics  The  system  characteristic  equation  for  the  y-channel  augmented  state 
equation  is  shown  in  Eq  (5.35).  A  Routh  stability  analysis  is  performed  to  determine 
stability(3)[185-191].  The  Routhian  array  for  the  characteristic  equation  is  shown  in  Table 
5.1. 


det  {SI  -  Ay)  =  +  (kyp  cos  a  -|- 1^  (kyp  +  kyi  cos  s  -f  fcyj  =  0  (5.35) 


(kyp  +  kyi  cos  a) 


T^w 

(kyp  COS  a  -I- 1)  kyi 

(fcypC08  0+l) 


"■yi 

Table  5.1.  Y-Channel  Routhian  Array 

After  applying  the  Routh  criterion,  the  following  stability  conditions  relating  the 
system  parameters  and  PI  controller  gains  are  obtained: 
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> 

0 

(5.36) 

k 

"’yp 

> 

1 

(5.37) 

cos  a 

kyp  -t-  kyi  cos  a 

> 

f^w  kyi 

(5.38) 

kyp  cos  a  -t- 1 

kyi 

> 

0 

(5.39) 

The  above  conditions  reveal  the  close  relationship  between  the  formation  geometry, 
controller  gains  and  stability  parameters.  This  relationship  is  critical  to  the  assignment 
of  the  PI  controller  gains.  In  the  next  section,  a  similar  analysis  is  performed  on  the 
x-channel. 

5.1.2  X-Channel  The  x-channel  is  decoupled  from  the  y-channel;  see  Eq  (5.23). 
The  decoupled  x-channel  state  equation  is  shown  below.  The  d{t)  signal  in  Equation  (5.40) 
is  considered  a  disturbance  input.  This  disturbance  perturbs  the  x-channel  due  to  a  wing 
heading  change  and  a  lead  aircraft  velocity  change;  the  disturbance  signal  is  defined  in 
Equation  (5.41). 


(5.40) 

d{t)  =  Vi  +  - ii’We  —  V’w)  (^-41) 

The  wing  aircraft’s  velocity  is  commanded  in  response  to  a  formation  separation 
change  in  the  x  direction.  The  PI  control  law  is  given  by 


where 


^Wc  —  -f- 


i: 


xdt 


(5.42) 
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kgp  =  non-dimensional,  x-separation  error,  proportional  gain 
kti  =  non-dimensional,  x-separation  error,  integral  gain 

Differentiating  Eq  (5.42)  yields 

Vwc  =  KiX  + 

Inserting  into  Eq  (5.43)  the  x  expression  from  Eq  (5.40)  gives 

V =  kgiX  kjgpV}^  kxpd 


(5.43) 


(5.44) 


The  augmented  state  equations  comprised  of  Eq  (5.40)  and  the  PI  controUer  (Eq  (5.44)), 
are  displayed  in  Eqs.  (5.45)  through  (5.48). 


i.e. 


1 


where 


Xx  =  AxXx  +  TxDx 

(5.45) 

Yx  =  CxXx 

(5.46) 

X 

0 

-1 

0 

X 

1 

K 

= 

0 

_ i_ 

fVyy 

1 

fvjy 

Vw 

+ 

0 

• 

kxi 

0 

Vwc 

kxp 

Y,= 


0  1  0 


(5.47) 


(5.48) 


Xj.  =  augmented  x-channel  state  vector 
Ax  =  augmented  x-channel  state  transition  matrix 
r,  =  augmented  x-channel  disturbance  input  matrix 
Dx  =  augmented  x-channel  disturbance  input  vector 
Yx  =  augmented  x-channel  output  vector 
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C,  =  augmented  x-channel  output  matrix 


Note:  The  primary  x-channel  parameter  in  Eq  (5.47)  is  the  velocity  time  constant  fv^,. 

(i)  Statics  Consider  a  unit  step  disturbance  input.  Similar  to  the  y-channel  “static” 
analysis,  the  steady  state  response  to  a  unit  step  disturbance  is 


f(t)„  =  lims(5*- A*)  =  = 

*-•0  s 


(5.49) 


Again,  the  wing  aircraft  must  be  able  to  track  the  lead  aircraft  with  zero  steady  state 
error.  This  zero  steady  state  tracking  error  requirement  is  stated  by 


At  i  —*  OQ,  i(oo)  =  0,  V(v(oo)  =  VWe(oo)  =  Vi  =  Dj,  =  1,  (5.50) 


i.e.. 


(5.51) 


Equation  (5.49)  implies  that  Eq  (5.51)  satisfies  this  requirement,  as  long  as  A*  is  invertible. 
Hence,  a  necessary  condition  for  the  augmented  state  matrix.  A,,  to  be  invertible  is  that 
the  integral  gain  must  be  non-zero.  This  is  in  line  with  the  results  from  the  y-channel 
analysis. 


(ii)  Dynamics  The  system  characteristic  equation  for  the  x-channel  augmented  state 
equation  is  shown  in  Eq  (5.52).  A  Routh  stability  analysis  is  again  performed  to  determine 
stability(.3)[l85-191].  The  Routhian  array  for  the  characteristic  equation  is  shown  in  Table 
5.2. 


det  (5/  -  A,)  =  fv^s^  +  +  k^pS  -|-  ib„-  =  0  (5.52) 

After  applying  the  Routh  criterion,  the  following  stability  conditions  on  both  the 
system  parameters  and  PI  controller  gains  are  obtained:  s 


5-11 


s 

s 


3 

2 


Table  5.2. 


TVyy  kxp 

1  ki 


X- Channel  Routhian  Array 


^Vw 

> 

0 

(5.53) 

> 

^xi'^Vw 

(5.54) 

kxi 

> 

0 

(5.55) 

The  above  equations  tie  in  the  type  of  aircraft,  controller  gains,  and  stability  param¬ 
eters.  This  relationship  is  critical  to  the  assignment  of  the  PI  controller  gains.  Hence  the 
following  holds: 

Proposition  1  The  necessary  and  sufficient  conditions  for  formation  flight  control  are  the 
employment  of  integral  feedback,  and  satisfaction  of  the  conditions  stated  in  Eqs  (5.36)- 
(5.39)  and  Eqs  (5.53)-(5.55).  □ 

The  next  section  makes  use  of  both  the  Y  and  X-Channels  analyses,  i.e.,  Eqs  (5.36)- 
(5.39)  and  Eqs  (5.53)-(5.55),  to  assign  the  optimal  PI  controller  parameters. 

5.2  Pole  Placement  through  Controller  Gain  Adjustment:  Residues  of  the 
Step  Response  for  Preliminary  Identification  of  Dominant  Dynamics 

The  closed-loop  MIMO  state  equations  developed  for  the  formation  flight  control 
problem  in  the  previous  section  are  used  to  design  the  PI  controller.  The  controller  is 
designed  in  two  steps.  In  the  first  step,  the  controller  gains  are  determined  for  the  linear 
state  equations  using  pole  placement.  The  second  step  involves  applying  these  gains  to  the 
non-linear  formation  control  system  model.  Computer  simulations  are  performed  to  fine 
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tune  these  gains  for  the  best  tracking  response.  The  poles,  residues  and  controller  gains 
presented  in  this  section  axe  all  dimensional. 

A  SISO  closed-loop  transfer  function,  Eq  (5.56),  is  obtained  for  both  the  Y  and  X- 
channels  by  optimizing  the  wing  aircraft’s  heading  and  velocity  response  respectively.  Using 
the  stability  and  steady  state  conditions,  from  the  previous  section,  a  transient  analysis 
is  performed  in  order  to  optimize  the  output  time  response.  The  dominant  dynamics  are 
identified  by  checking  the  residues  of  the  output  response.  The  controller  gains  are  then 
adjusted  to  produce  the  fastest  poles  with  the  largest  residues. 


G(s)  =  Ce,(s/-A„]-'r,,  (5.56) 

Computer  simulations  using  the  non-linear  control  system  are  performed  using  these 
gains.  The  rate-limiters  in  the  aircraft  models  tend  to  produce  an  underdamped  oscillatory 
response  for  high  integral  gains  values.  This  response  is  caused  by  the  increase  in  tracking 
error  induced  by  the  presence  of  the  rate-limiters.  Fiqure  5.2  shows  an  example  of  this  build 
up.  Indeed,  integral  action  is  needed  for  tracking.  At  the  same  time,  integral  feedback 
adversely  interacts  with  the  saturation.  Hence,  an  interesting  trade  off  situation  has  arisen, 
which,  however,  is  not  fully  investigated  in  this  thesis. 

In  the  next  two  sections,  this  analysis  is  conducted  on  both  the  Y  and  X  channels 
in  order  to  develop  an  optimum  set  of  gains.  This  analysis  is  performed  using  the  C-130A 
first  order  models  developed  by  Dargan(2)  and  Rohs(9).  The  two  aircraft  are  arranged  in 
a  diamond  formation  as  shown  in  Figure  5.1.  The  key  dimensional  parameters  needed  for 
this  analysis  are  contained  in  Table  5.3. 
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Table  5.3.  Formation  and  Aircraft  Parameters 


Parameter 

Description 

Value 

Units 

V 

Nominal  Velocity 

375 

fps 

Xo 

x-separation 

500 

feet 

Vo 

y-separation 

500 

feet 

a 

formation  angle 

45 

degs 

1 

separation 

707 

feet 

Hdg  time  constant 

1 

}  ® 

secs 

TVw 

Vel  time  constant 

1 

_a _ 

secs 
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5.2.1  Y-Channel  Gains  The  Y-Channel  controller  gains  are  selected  by  optimiz¬ 
ing  the  wing  aircraft’s  heading  response  to  a  lead  aircraft’s  step  heading  change.  The 
Laplace  representation  of  the  output  response  is 


*„W  =  <3WfeW=[o  1  o][»/-/l,l-T,Ul  +  ^  +  ^  +  ^  (5.57) 


Both  kyi  and  kyp  are  varied,  and  the  resulting  poles  and  residues  of  Eq  (5.57)  are  determined 
using  computer  simulations.  The  ratio  of  is  fixed,  while  kyp  is  varied  from  0.001  to  10. 
Figure  5.3  displays  the  poles  and  residues  for  =  100.  The  best  response  is  attained  for 
values  of  kyp  >  0.1.  The  pole  is  no  longer  the  dominant  pole  since  its  associated  residue 
is  nearly  zero. 


.001  .01  .1  1  10 

•«yp 


Figure  5.3.  Y-Channel,  Wing’s  Heading  Step  Response,  Poles  and  Residues  versus  Con¬ 
troller  Gains,  ^  =  100 


Selecting  kyp  =  1  and  kyi  =  0.01  results  in  the  following  poles  and  residues: 
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.0010 


(5.58) 


,  ,  ,  1  ,  0.0001  ,  -1.0010 
V’iv(s)  -  -  +  ^  ^  ^  ^  .7492  5  +  750.7408 


The  output  time  response  is  obtained  through  the  inverse  Laplace  transform  of 
Eq  (5.58).  Other  ratios  of  may  be  selected  but  numerous  computer  simulations  have 
shown  this  to  be  a  near  optimum  value. 


5.2.2  X-Channel  Gains  The  x-channel  controller  gains  are  selected  by  optimizing 
the  wing  aircraft’s  response  to  a  lead  aircraft’s  step  velocity  change.  The  disturbance 
input  in  Eq  (5.41)  is  confined  to  the  lead  aircraft’s  velocity  by  setting  both  the  heading 
and  commanded  heading  to  zero.  The  wing’s  velocity  response  to  a  step  input  is  shown 
inEq  (5.59). 


V.(,)  =  G.(»)C,(»)  =  [010]  [»/  -  =  1  +  ^  ^  ^  (5.59) 

The  poles  and  residues  of  Eq  (5.59)  are  again  determined  using  computer  simulation. 
The  ratio  of  ^  is  fixed  while  k^p  is  varied  from  0.001  to  10.  Figure  5.4  displays  the  poles 
and  residues  for  =  100.  A  fast  response  is  attained  for  values  of  k^p  >  0.7.  The  pole 
P2  is  not  dominant  since  it’s  associated  residue  is  much  less  than  the  other  three  residues. 

Selecting  k^p  =  0.75  and  fc*,-  =  0.0075  results  in  the  following  output  poles  and 
residues: 


,  1  0.0138  -6.6677  .  5.6538 

Viy(s)  —  —  “i“ - f“  -  "1“  - 

The  next  section  compares  the  controller  gains  determined  using  the  above  analysis 
against  those  developed  in  previous  research(2)[4- 1-4-13].  This  comparison  is  done  using 
the  non-linear  formation  control  system  model. 
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•w 


Figure  5.4.  X-Channel,  Wing’s  Velocity  Step  Response,  Poles  and  Residues  versus  Con¬ 
troller  Gains,  ^  =  100 

Kyi 


5.3  Performance  Evaluations 

Two  different  computer  simulations  are  performed  to  evaluate  the  controller  designed 
in  the  previous  two  sections.  A  diamond  formation  consisting  of  a  single  wing  aircraft 
is  commanded  to  have  first  a  velocity  change,  and  then  a  heading  change.  The  wing 
aircraft’s  time  response  is  compared  against  values  developed  in  previous  research(2)[4-l- 
4-13].  Table  5.4  summarizes  the  PI  controller  values  used  in  the  evaluations.  No  mixer  is 
used  in  this  evaluation. 


Table  5.4.  PI  Controller  Gains 


Parameter 

New 

Previous 

kyp 

1 

.2 

kyi 

.01 

.009 

kgp 

.75 

.2 

kxi 

.0075 

.015 
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5.3.1  Diamond  Formation  Heading  Change  The  response  for  a  10  degree 
heading  change  from  a  diamond  formation  is  shown  in  Figure  5.5.  This  figure  shows  a 
slightly  ftister  response  of  the  wing  aircraft  in  tracking  the  lead  aircraft.  There  is  a  sig¬ 
nificantly  reduced  error  in  the  y-separation  during  the  maneuver.  The  better  tracking  is 
offset  by  a  slight  increase  in  the  wing’s  velocity. 


rtow  -  aolid,  PravkMia  -  dottad 


Figure  5.5.  Diamond  Formation,  10  Degrees  Heading  Change,  Time  Response 
Comparison 


5.3.2  Diamond  Formation  Velocity  Change  The  response  for  a  400  feet  per 
second  velocity  command  is  shown  in  Figure  5.6.  This  figure  shows  a  faster  response  of  the 
wing  aircraft  in  tracking  the  lead  aircraft.  The  new  controller  gains  result  in  virtually  zero 
overshoot,  since  the  output  poles  are  overdamped  and  the  integral  gain  is  relatively  small. 
This  low  integral  gain  adversely  affects  the  y-separation,  resulting  in  a  high  settling  time. 
In  terms  of  the  velocity  response,  the  new  gains  definitely  result  in  a  better  response.  The 
integral  gain  k^i  could  be  increased  to  decrease  the  y-separation  settling  time. 
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N«w  -  solid,  Previous  -  dottsd 


Figure  5.6.  Diamond  Formation,  375  to  400  fps  Velocity  Change,  Time  Response 
Comparison 


5.4  Summary 

An  analytical  analysis  of  the  formation  flight  control  problem  produced  a  linear  four- 
dimensional  state  space  model  using  small  angle  approximations  and  small  perturbations. 
The  X  and  Y-Channels  are  decoupled  and  augmented  using  a  PI  control  law.  Further 
analysis  of  these  augmented  closed-loop  equations  substantiated  the  requirement  for  inte¬ 
gral  feedback  to  achieve  zero  steady-state  error.  Application  of  Routh’s  stability  criterion 
identified  the  interrelationship  of  formation  geometry,  aircraft  models,  and  controller  gains 
on  stability.  Simple  pole  placement  is  used  to  select  an  optimum  set  of  PI  controller  gains. 
The  performance  evaluation  showes  that  these  gains  produce  a  fast  time  response  with  zero 
steady  state  error.  While  the  analytical  analysis  result  in  acceptable  PI  controller  gains, 
additional  refinement  of  these  values  is  still  needed  to  produce  the  best  tracking  response 
in  the  non-linear  computer  simulation. 


5-19 


VI.  Energy  Conserving  Maneuvers 


Horizontal  and  constant  velocity  lead  aircraft  maneuvers  force  the  wing  aircraft  to 
adjust  its  speed  to  maintain  the  formation.  Consequently,  the  wing  aircraft  is  constantly 
accelerating  and  decelerating  during  level  fUght  formation  turns.  Conserving  the  wing 
aircraft’s  energy,  will  reduce  fuel  consumption.  In  addition,  it  is  advantageous  from  a 
reliability  and  maintainability  point  of  view  to  have  the  engine  operate  at  a  constant 
throttle  setting.  Hence,  the  vertical  separation  requirement  between  the  lead  and  the  wing 
aircraft  is  relaxed  in  order  to  allow  for  energy  conserving  maneuvers.  This  allows  the  wing 
aircraft  to  dive  to  increase  its  velocity,  or  to  zoom  to  slow  down.  The  next  section  presents 
the  theoretical  analysis  outlining  a  wing  energy  conserving  formation  control  strategy.  As 
in  the  previous  chapter,  first  order  aircraft  models  are  used.  A  performance  evaluation  of 
the  energy-conserving  formation  flight  control  system  follows  this  analysis. 

6.1  Theoretical  Development 

The  wing  aircraft’s  specific  energy  is  the  sum  of  its  kinetic  and  potential  energy  per 
unit  mass: 

cw  =  2^w  +  (6-1) 

where 

g  =  acceleration  of  gravity 

ew  =  wing  aircraft’s  specific  energy 

Diflferentiating  Eq  (6.1)  and  setting  eyy  equal  to  zero  results  in  Eq  (6.2).  This  equation  is 
the  fundamental  relationship  for  energy  conservation. 


evv  =  VVw  -b  ghw  =  0 


(6.2) 


Rearranging  this  equation  and  substituting  the  first  order  aircraft  model 


:  1  .  1  . 
hw  — - nw  H - hwc 


(6.3) 
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yields 


(6.4) 


hwc 


As  seen  in  Eq  (6.4),  the  wing  aircraft’s  velocity  is  controlled  by  using  an  altitude  change 
command. 


The  following  equations  display  the  definitions  needed  to  non-dimensionalize  Eq  (6.4). 


hw  —  ^  hwe  = 


Furthermore, 


Vw  —  ir—hw  —  -:^hwc 


and  the  nondimensional  parameter 


ii 


where 


(6.5) 


(6.6) 


(6.7) 


^  =  altitude  to  velocity  proportionality  constant 


A  PI  controller  is  designed  to 
x-separation.  The  appropriate 
defined  in  Eq  (5.41). 


command  the  wing’s  altitude  in  response  to  a  change  in 
equations  are  listed  below.  The  disturbance  input  d{t)  is 


‘Wc 


=  KzpX  +  kx 


i: 


xdt 


(6.8) 


where 


kxzp  =  non-dimensional,  3-D  x-separation  error,  proportional  gain 
kxzi  =  non-dimensional,  3-D  x-separation  error,  integral  gain 

Differentiating  Eq  (6.8)  yields: 


=  k^^iX  -1-  k 


(6.9) 
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Inserting  the  x  from  Eq  (5.6)  and  d{t)  from  Eq  (5.42)  into  Eq  (6.8)  yields 


=  kcziX  +  kxjip^hw  +  kxzpd 


(6.10) 


The  dosed  loop  X  channel  state  space  equations  incorporating  energy  conservation  and 
proportional  plus  integral  feedback  are  shown  in  Eqs  (6.11)  to  (6.12). 


where 


X  —  AxzXxz  + 


(6.11) 
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(6.12) 


(6.13) 


Xxz  =  augmented  x-channel  state  vector 
Axz  =  augmented  x-channel  state  transition  matrix 
=  augmented  x-channel  disturbance  input  matrix 
Dxz  =  augmented  x-channel  disturbance  input  vector 
Yxz  =  augmented  x-channel  output  vector 
Cxz  =  augmented  x-channel  output  matrix 


Given  a  unit  step  disturbance  input,  the  wing  aircraft  must  be  able  to  track  the  lead  with 
zero  steady  state  error.  The  steady  state  response  to  a  unit  step  disturbance  is 


lims(5/->l„)-‘r„i 


— >4“^r  — 

^xz  * 


(6.14) 
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Zero  steady  state  error  is  verified  by  rearranging  Eq  (6.14).  Equation  (6.15)  verifies  that 
this  requirement  is  satisfied  as  long  as  is  invertible. 


0 

-1 

r  —A 

^  XZ  ^xz 

1 

= 

0 

1 

-kzi 

(6.15) 


Analysis  of  the  matrix  in  Eq  (6.12)  shows  that  a  necessary  condition  for  the  augmented 
state  matrix,  A^z,  to  be  invertible  is  that  the  integral  gain  must  be  non-zero.  This 
reaffirms  Proposition  1,  in  Chapter  V,  that  integral  feedbaurk  is  a  necessary  and  sufficient 
condition  for  formation  ffight  control. 

The  system  characteristic  equation  for  the  X-Channel  augmented  state  equation  is 
shown  is  Eq  (6.16).  A  Routh  stability  analysis  is  performed  to  determine  stability.  The 
Routhian  array  for  the  characteristic  equation  is  shown  in  Figure  6.1. 


det  {si  -Az.z)  =  +  5^  -  Kzpis  -  kz,zi^  =  0  (6.16) 

After  applying  the  Routh  criterion,  the  following  stability  constraints  on  both  the  system 


kxzp^ 

1 

^~kzp^  + 

kzi^ 


Table  6.1.  X-Channel  Routhian  Array 


parameters  and  PI  controller  gains  result. 


kzi 


>  0 
k 


kzi  <  0 


(6.17) 

(6.18) 
(6.19) 


The  above  equations  impose  the  requirement  for  negative  controller  gains  in  order  to 
maintain  stability.  This  inverse  relationship  is  easily  visualized,  since  the  wing  aircraft 
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must  descend  to  speed  up  and  climb  to  slow  down.  The  next  section  uses  an  analysis  of 
the  X-Channel  to  assign  the  optimum  PI  controller  parameters. 


6.2  X-Channel  Controller  Gain  Adjustment 

The  X-Channel  controller  gains  are  selected  by  optimizing  the  wing  aircraft’s  response 
to  a  lead  aircraft’s  velocity  change.  When  the  heading  and  the  commanded  heading  are 
equal,  the  disturbance  input  given  by  Eq  (5.41)  is  reduced  to  the  lead  aircraft’s  velocity. 
The  wing’s  velocity  response  to  a  step  input  is 


^w{s)  =  0  - 


i  0 


[s/  -  A„]  —  -  + 


+ 


’*3 


+ 


Ta 


S  +  P2  s+ps  s+ Pa 


(6.20) 


The  poles  and  residues  of  Eq  (6.20)  are  again  determined  using  computer  simulation.  The 
ratio  of  is  fixed  while  k^p  is  varied  from  0.001  to  100.  Figure  6.1  displays  the  poles  and 
residues  for  =  100.  A  fast  response  is  attained  for  a  value  of  k^ip  <  —2.  Again,  the 
pole,  p2,  is  a  less  dominant  pole  since  it’s  associated  residue  is  much  smaller  less  than  the 
residues  of  the  other  poles.  Selecting  =  —3  and  kx,i  =  -0.03  results  in  the  response 
shown  in  Eq  (6.21).  Even  though  the  selected  poles  are  complex,  they  provide  a  faster 
response  than  selecting  real  poles. 


1  0.0432  1.4561/1-  135.7624  deg  1.4561/1+  135.7624  deg  , 

-  + - - - ^ - -  (6.21) 

s  s  +  0.0104  s  +  0.2448  -  j0.2510  s  +  0.2448  +  j0.2510  ''  ’ 


In  the  next  section,  a  set  of  controller  gains  which  produces  a  good  response  is  determined, 
using  the  above  analysis,and  through  computer  simulation. 


6.3  Performance  Evaluations 

Two  different  computer  simulations  are  performed  to  evaluate  the  controller  designed 
in  the  previous  section.  A  diamond  formation  consisting  of  a  single  wing  aircraft  is  com¬ 
manded  to  perform  a  a  30  degree  “side  step”  heading  change.  The  second  maneuver  is  a 
90  degree  heading  and  formation  change.  Table  6.2  summarizes  the  PI  controller  values 
used  in  the  evaluations.  No  mixer  is  used  in  this  evaluation. 
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Figure  6.1.  X-Channel  Design  for  Energy  Conservation,  Wing’s  Velocity  Step  Response, 
Poles  and  Residues  versus  Controller  Gains,  =  100 


6.3.1  Diamond  Formation  Side  Step  Maneuver  The  side  step  heading  change 
maneuver  consists  of  two  consecutive  equal  but  opposite  lOdeg  heading  commands.  The 
purpose  of  this  maneuver  is  to  return  the  formation  back  to  its  original  heading  after  the 
initial  command.  Figure  6.2  displays  the  time  responses  for  this  maneuver.  During  the 
initial  lOdeg  heading  change,  the  wing  aircraft  must  descend,  increasing  its  velocity,  in 
order  to  track  the  lead.  As  the  lead  aircraft  executes  the  -lOdeg  heading  change  during 
the  second  part  of  the  maneuver,  the  wing  ascends,  bleeding  off  its  airspeed. 


Table  6.2.  PI  Controller  Gains 


Parameter 

Value 

kyp 

1 

kyi 

0.01 

-0.3 

krii 

-0.03 
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Figure  6.2.  Diamond  Formation  10  Degree  Side  Step  Heading  Change 


6.3.2  Diamond  Formation  90  Degree  Heading  Change  and  Formation 
Transposition  As  seen  in  the  side  step  heading  change,  the  wing  aircraft  must  increase 
its  velocity  in  order  to  follow  the  lead  aircraft  for  a  positive  heading  change.  If  the  lead 
is  to  make  a  90  degree  heading  change,  a  considerable  amount  of  energy  is  required  to 
increase  the  wing’s  velocity  and  track  the  lead  aircraft  along  the  outside  path.  Fiqure  6.3 
illustrates  an  alternative  trajectory  for  the  wing  aircraft.  In  this  case,  the  wing  aircraft 
conserves  its  energy  by  decreasing  its  altitude,  which  increases  it’s  velocity,  and  then  exe¬ 
cutes  a  formation  change  to  the  opposite  side  of  the  lead  aircraft  as  is  often  an  operational 
practice.  The  wing  then  climbs  to  bleed  off  the  excess  velocity  and  resumes  formation 
holding  on  the  other  side  of  the  lead. 

A  design  is  accomplished  in  which  the  gains  are  changed  from  those  in  Table  7.1. 
The  PI  controller  gains  are  reduced  to  kg^p  =  —0.5  and  kg,i  =  -0.005  to  facilitate  this 
maneuver.  The  time  responses  for  this  maneuver  are  shown  in  Figure  6.4.  Figure  6.5 
displays  the  inertial  flight  paths  for  both  aircraft.  The  circular  dots  are  placed  at  five 
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Figure  6.3.  Diamond  Formation  90  Degree  Heading  and  Formation  Change  utilizing  En¬ 
ergy  Conservation,  Flight  Path  Hlustration 


second  intervals.  Even  though  the  wing’s  path  crosses  the  lead’s,  there  is  no  collision 
because  they  are  at  different  altitudes.  This  is  illustrated  in  the  final  plot,  which  displays 
a  three  dimensional  view  of  this  maneuver. 
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Figure  6.4.  Diamond  Formation  90  Degree  Heading  and  Formation  Change  Utilizing  En¬ 
ergy  Conservation;  Time  Responses 
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Figure  6.5.  Diamond  Formation  90  Degree  Heading  and  Formation  Change  Utilizing  En¬ 
ergy  Conservation,  Flight  Path  in  Inertial  Frame 
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Figure  6.6.  Diamond  Formation  90  Degree  Heading  and  Formation  Change  Utilizing  En¬ 
ergy  Conservation;  3-D  Flight  Path  in  Inertial  Frame 

6. 4  Summary 

This  chapter  introduces  the  concept  of  en>?rgy  conservation  into  formation  flight  con¬ 
trol.  By  maintaining  the  wing  aircraft’s  energy  constant,  fuel  consumption  is  reduced 
during  horizontal  maneuvers.  Linearized  equations  are  developed  that  facilitate  the  syn¬ 
thesis  of  the  control  law  needed  to  implement  this  concept.  A  PI  controller  is  used  to 
command  the  wing  aircraft  in  response  to  separation  changes  caused  by  the  lead  aircraft. 
The  controller  successfully  maintains  the  formation  through  a  side  step  maneuver  while 
maintaining  the  wing’s  energy  constant.  A  second  maneuver  consisting  of  a  90  degree 
heading  change  is  also  simulated.  An  alternate  trajectory  is  chosen  for  the  wing  aircraft 
in  order  to  conserve  energy. 
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VII.  Automated  Formation  Flight  Control  Simulation  using 
Second-Order  C-130  Aircraft  Models 

In  this  chapter,  the  second  order  aircraft  models  developed  in  Appendix  B  are  used 
in  the  computer  simulation.  The  formation  flight  control  system  as  described  in  Chapter 
II,  Section  2.5.2,  implementing  a  PI  controller  and  linear  mixer  is  evaluated.  This  control 
system  was  designed  in  previous  research  using  first  order  aircraft  models.  The  purpose  of 
this  chapter  is  to  evaluate  the  effects  on  the  system  response  due  to  the  new  aircraft  models. 
Of  special  interest  is  the  transient  time  response  of  the  wing  aircraft  during  a  formation 
maneuver.  The  new  aircraft  models  incorporate  an  onset  delay  in  both  the  heading  and 
the  altitude  channels.  Both  the  PI  controller  and  linear  mixer  gains  need  to  be  adjusted  in 
order  to  compensate  for  this  time  delay.  In  the  next  section,  the  formation  control  system  is 
simulated  using  the  PI  controller  and  gains  developed  in  Dargan’s  research  (2)  [4- 14-4-22]. 

7.1  Initial  Performance  Evaluation 

The  formation  flight  control  system  utilizing  the  PI  controller  and  linear  mixer,  as 
illustrated  in  Fiqure  2.7,  produce  a  fast  time  response  and  zero  steady-state  error,  when 
simulated  using  first  order  aircraft  models(2)[Fig  4.12].  The  new  aircraft  models,  as  shown 
in  Figure  2.2,  display  a  slower  time  response  than  the  first  order  models.  Even  though  the 
velocity  channel  remains  a  first  order  transfer  function,  its  time  constant  is  much  laxger 
than  in  the  previous  model.  This  results  in  a  longer  time  response.  Computer  simulations 
consisting  of  formation  velocity,  heading  and  altitude  change  are  evaluated.  Table  7.1  lists 
the  key  parameters  used  in  the  following  simulations.  Note,  the  mixer  is  not  used  for  the 
altitude  channel. 

7.1.1  Velocity  Change  Performance  Evaluation  The  simulation  as  described 
in  Table  7.1  is  commanded  a  step  velocity  change  input  from  a  nominal  velocity  of  375 
feet  per  second  to  400  feet  per  second.  The  time  response  plots  for  this  input  are  shown 
in  Figure  7.1.  While  the  lead  aircraft’s  velocity  response  exhibits  a  first  order  overdamped 
response,  the  wing’s  response  is  clearly  second  order  and  underdamped.  The  x-separation 
is  also  underdamped,  and  has  a  large  overshoot.  Clearly  these  responses  are  less  desirable 
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Table  7.1.  Initial  Performance  Evaluation  Parameters 


Parameter 


Lead  Aircraft 


Wing  Aircraft 


PI  Controller  Gains 


X  Channel 


Y  Channel 


Z  Channel 


Linear  Mixer  Gains 


X  Channel 


Y  Channel 


Z  Channel 


Formation 


X  separation 


y  separation 


y  separation 


Nominal  Velocity 
hline  Nominal  Heading 


Nominal  Altitude 


Value 


C-130B 


C-130B 


Diamond 


500  feet 


500  feet 


0  feet 


375  fps 
0  degrees 


500  feet 


than  the  time  responses  using  first  order  aircraft  models.  On  the  positive  side,  the  system 
is  stable,  with  zero  steady  state  error.  The  previous  aircraft  model  time  responses  are 
shown  as  dotted  lines. 


7.1.2  Heading  Change  Performance  Evaluation  The  formation  control  sys¬ 
tem  response  to  a  30  degree  heading  change  is  displayed  in  Figure  7.2.  Again  the  second 
order  aircraft  models  result  in  a  less  than  optimal  time  responses.  Even  though  there  is 
minimal  degradation  in  the  wing’s  heading  response,  the  wing’s  velocity  and  the  x  and 
y  separation  responses  are  highly  oscillatory.  Again,  on  the  positive  side,  the  system  is 
stable  with  zero  steady  state  error. 
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Time  Response  Plots  to  a  400  fps  Velocity  Command,  1st  and  2nd  Order 
Aircraft,  Initial  Gains 


Figure  7.2.  Time  Response  Plots  to  a  30  Degrees  Heading  Command,  1st  and  2nd  Order 
Aircraft,  Initial  Gains 
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7.1.3  Altitude  Change  Performance  Evaluation  The  formation  control  sys¬ 
tem  response  to  550  feet  altitude  change  is  shown  in  Figure  7.3.  The  altitude  response 
is  clearly  the  best  of  the  three  channels.  There  is  a  slightly  noticeable  oscillation  present 
in  both  the  wing’s  altitude  response  and  z  separation.  This  response  compares  favorably 
with  the  first  order  system  response.  Consequently,  the  system  is  stable  with  zero  steady 
state  error. 
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Figure  7.3.  Time  Response  Plots  for  a  550  feet  Altitude  Change  Command,  1st  and  2nd 
Order  Aircraft,  Initial  Gains 
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1.2  Performance  Evaluation  using  Modified  PI  Controller  and  Linear  Mixer 

Gains 

The  computer  simulations  in  the  previous  section  clearly  indicate  that  the  new,  pre¬ 
dominantly  second-order  and  slower  aircraft  models,  yield  formation  control  system  re¬ 
sponses  that  are  not  as  good  as  those  provided  by  the  first-order  aircraft  models.  The 
mixer  and  controller  parameters  determined  using  the  first  order  aircraft  models  no  longer 
produce  satisfactory  performance  using  the  slower,  second  order  models.  Consequently, 
the  PI  controller  and  mixer  gains  are  adjusted  to  improve  the  formation  control  system 
time  response.  The  controller  gains  are  adjusted  in  order  to  reduce  the  oscillations  and 
instability  caused  by  the  previous  values.  These  new  values  are  shown  in  Table  7.2. 


Table  7.2.  Adjusted  PI  Controller  and  Mixer  Performance  Parameters 


Parameter 

Value 

PI  Controller  Gains 

X  Channel 

Y  Channel 

Z  Channel 

Linear  Mixer  Gains 

X  Channel 

A:,  =  1,  =  5 

Y  Channel 

=  1,A:^  =  0 

Z  Channel 

A:^  =  1 

Both  the  X  and  Y-Channel  gains  are  increased  significantly  in  order  to  optimize 
the  response.  The  linear  mixer  is  the  key  in  smoothing  out  the  wing’s  velocity  response. 
However,  the  heading  error  gain,  A:^,  is  disabled  since  it  has  very  little  effect  on  the  time 
response.  The  integral  PI  controller  gain  in  the  Z  channel  is  reduced  in  order  to  dampen  the 
response.  The  same  computer  simulations  performed  in  the  previous  section  are  repeated 
using  these  new  parameters. 

1.2.1  Velocity  Change  Performance  Evaluation  The  nominal  formation  pa¬ 
rameters  as  described  in  Table  7.1,  in  addition  to  the  PI  controller  and  linear  mixer  gains 
listed  in  Table  7.2,  are  used  in  this  evaluation.  Again  a  step  velocity  change  input  from 
a  nominal  velocity  of  375  feet  per  second  to  400  feet  per  second  is  simulated.  The  time 
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response  plots  for  this  input  are  shown  in  Figure  7.4.  The  previous  first  order  responses 
are  again  repeated  on  the  same  plot  for  comparison.  Even  though  the  new  aircraft  model 
responses  are  naturally  slower  than  the  previous  models,  they  exhibit  smooth  response 
with  no  overshoot.  As  expected,  the  new  aircraft  models  do  not  reach  the  rate  limiters  as 
is  clearly  the  case  in  the  old  models. 


Figure  7.4. 


Time  Response  plots  of  a  400  fps  Velocity  Command,  Diamond  Formation, 
1st  and  Second  Order  Models  and  Initial  Controller  Gains 


1.2.2  Heading  Change  Performance  Evaluation  The  formation  response  to 
a  30  degree  heading  change  is  shown  if  Figure  7.5.  After  implementing  the  new  set  of  PI 
controller  gains,  both  the  first  and  second  order  system  responses  are  very  similar.  Both 
the  lead  and  wing  aircraft  saturate  their  respective  rate  limiters  during  the  maneuver. 
Consequently,  their  time  responses  are  similar.  This  may  indicate  an  inconsistency  between 
the  rate  limit  values,  and  tli'’i  transfer  functions  used  to  model  the  heading  channel.  The 
y-separation  response  for  the  second  order  model  displays  a  high  frequency  oscillation  not 
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seen  in  the  heading  response.  Changing  the  proportional  or  integral  gains  is  effective  in 
dampening  this  oscillation. 
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Figure  7.5.  Time  Response  plots  of  a  30  Degree  Heading  Command,  Diamond  Formation, 
1st  and  Second  Order  Models  and  Modified  Controller  Gains 


7.2.3  Altitude  Change  Performance  Evaluation  The  formation  control  sys¬ 
tem  response  to  a  550  feet  altitude  change  command  is  shown  in  Figure  7.3.  The  slight 
decrease  in  k^i  removes  the  light  oscillation  present  in  the  earlier  second  order  response. 
Other  than  the  speed  of  the  response,  both  altitude  models  exhibit  similar  time  responses. 


7.3  Summary 

This  chapter  provides  a  preliminary  evaluation  of  the  performance  of  the  formation 
flight  control  system  utilizing  the  new,  second-order  and  slower  aircraft  models.  Computer 
simulation  of  three  separate  formation  maneuvers  illustrated  the  effects  of  the  new  and 
“slower”  aircraft  models  on  the  control  system  response.  Both  the  velocity  and  heading 
channels  display  an  oscillatory  response  with  poor  settling  time.  The  least  degradation 
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Figure  7.6. 


Time  Response  plots  of  a  550  feet  Altitude  Command,  Diamond  Formation, 
1st  and  Second-Order  Models  and  Initial  Controller  Gains 


is  in  the  altitude  channel.  Consequently,  the  same  control  system  parameters  can  not  be 
used  for  first  and  second-order  models.  An  axljustment  of  both  the  PI  controller  and  mixer 
gains  results  in  satisfactory  transient  responses  with  zero  steady  state  error.  Except  for 
the  heading  channel,  the  new  aircraft  models  exhibit  satisfactory  and  more  realistic  time 
responses. 
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VIIL  Analysis  and  Conclusions 


8.1  Analysis  of  Results 

The  objective  of  this  research  is  to  continue  the  development  of  an  automated  forma¬ 
tion  flight  control  system  through  an  analytical  analysis  of  the  control  problem,  investiga¬ 
tion  of  energy  preserving  formation  flight  control  strategies,  development  of  more  realistic 
aircraft  models,  improved  visualization  tools,  and  computer  simulation.  These  objectives 
have  been  met. 

8.2  Display  of  Lead  Aircraft’s  Trajectory  as  a  Means  of  Evaluating  the  For* 

mation  Flight  Control  System’s  Performance 

Computer  simulation  is  used  to  evaluate  the  design  of  the  formation  control  system. 
Normally,  the  time  response  plots  of  each  aircraft’s  velocity,  heading,  altitude,  and  sepa¬ 
rations  distances  are  evaluated  after  each  computer  simulation.  As  shown  in  Chapter  III, 
a  graphical  method  of  displaying  the  lead  jurcraft’s  trajectory  in  the  wing  aircraft’s  refer¬ 
ence  frame  is  developed.  Analogous  to  Lissajous  flgures  used  to  visualize  two  dimensional 
oscillations,  this  method  is  a  valuable  tool  in  evaluating  the  tightness  and  speed  of  the 
formation  control  system. 

8.3  Two  Wing  Aircraft  Formation  Flight  Control 

Using  first-order  aircraft  models,  an  aircraft  formation  consisting  of  two  wing  aircraft, 
is  simulated  in  Chapter  IV.  Each  wing  aircraft  uses  both  PI  control  and  a  mix  of  separation 
and  maneuver  errors  to  track  the  leader.  The  wing  aircraft  successfully  tracks  the  lead 
aircraft  through  30  degree  formation  heading  changes.  No  coordination/communication 
between  the  wing  aircraft  is  used.  Consequently,  multiple  wing  aircraft  formations  are 
feasible  as  long  as  the  formation  separation  distance  restrictions  are  maintained.  These 
restrictions  are  based  on  each  aircrafts’s  tracking  error  which  are  determined  using  single 
wing  aircraft  computer  simulations. 
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8.4  Analytical  Formation  Flight  Controller  Synthesis 

In  previous  research  efforts,  the  key  parameters  in  the  formation  control  system  were 
developed  using  a  trial  and  error  method.  These  parameters  included  the  PI  controller 
and  linear  mixer  gains.  This  method  produced  successful  results,  yet  required  considerable 
effort  and  numerous  computer  simulations. 

An  analytical  approach  to  determining  the  optimal  control  gains  is  developed  in 
Chapter  V.  Using  first-order  aircraft  models,  a  linearized  and  parameterized  MIMO  plant 
is  developed  and  augmented  using  the  PI  control  law.  Steady  state  analyses  of  the  closed- 
loop  state  matrix  reveals  the  necessary  and  sufficient  condition  of  integral  feedback  control 
for  formation  ffight  control.  A  Routh  stability  analysis  on  the  closed  loop  characteristic 
equation  shows  the  close  relationship  between  the  formation  geometry,  controller  gains  and 
system  stability. 

A  two  step  process  is  developed,  using  the  above  analysis,  to  determine  the  PI  con¬ 
troller  gains.  In  the  first  step,  pole  placement  is  used  to  select  a  preliminary  set  of  controller 
gains.  Computer  simulations  of  the  non-linear  formation  ffight  control  system  are  then 
performed  in  order  to  fine-tune  these  preliminary  gains.  Fine  tuning  through  computer 
simulation  is  still  needed  to  take  into  account  the  non-linear  rate  limiters  present  in  the 
aircraft  models.  When  compared  with  the  controller  gains  developed  using  a  purely  trial 
and  error  approach,  these  new  values  perform  somewhat  better.  The  development  of  the 
linear  closed  loop  state  space  equations  provides  valuable  insight  into  the  formation  ffight 
control  problem  and  aids  in  determination  of  PI  controller  gains. 

8.5  Energy  Conservation  of  Wing  Aircraft 

The  primary  need  for  energy  conservation  of  the  wing  aircraft  is  to  reduce  fuel  con¬ 
sumption  during  large  formation  heading  changes.  If  the  wing  is  positioned  in  a  diamond 
formation,  an  initial  increase  in  velocity  is  required  to  track  the  lead  along  the  “outer 
track”  trajectory  during  a  heading  change.  This  is  followed  by  a  speed  reduction  as  the 
wing  mrcraft  settles  back  into  formation  once  the  maneuver  is  complete.  The  opposite 
is  true  for  the  “inside  track”  trajectory  of  the  diamond  formation.  Instead  of  using  the 
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throttle  to  first  increase  and  then  decrease  velocity,  the  wing  aircraft  descends  to  increase 
velocity,  and  then  ascends  to  decrease  its  velocity,  while  preserving  its  total  energy. 

An  analytical  approach  is  used  to  design  a  formation  flight  control  system  using  an 
altitude  change  to  conserve  the  wing’s  energy  during  velocity  changes.  A  formation  control 
system  based  on  a  PI  controller  and  first-order  aircraft  models  is  designed.  Pole  placement 
through  analysis  of  the  systems  step  response  is  used  to  select  the  controller  gains.  Com¬ 
puter  simulations  show  that  the  control  system  maintains  the  formation  throughout  the 
heading  maneuver,  while  maintaining  the  wing’s  energy  constant.  Wing  transpositions 
during  heading  change  maneuvers  are  also  considered. 

8. 6  Computer  Simulation  of  Formation  Flight  Control  System  using  Second- 

Order  Aircraft  Models 

Up  to  this  point,  first-order  aircraft  autopilot  models  have  been  used  to  model  the 
flight  performance  of  the  C-130  aircraft.  The  development  of  second-order  aircraft  models 
is  addressed  in  Appendices  A  and  B.  New  aircraft  models  based  on  computer  simulation  of 
C-130  longitudinal  and  lateral  autopilots  are  designed.  The  new  models  incorporate  second- 
order  transfer  functions  in  both  the  velocity  and  heading  channels.  The  second-order 
transfer  functions  better  model  the  time  delay  present  in  the  autopilot’s  step  response.  In 
addition,  the  new  and  “slower”  aircraft  models  exhibit  a  more  realistic  and  linear  behavior, 
by  not  immediately  saturating  the  rate  limits  during  command  inputs. 

In  Chapter  VII,  computer  simulation  of  velocity,  heading  and  altitude  formation 
maneuvers  show  that  the  new  models  affect  the  response  of  the  control  system  as  expected. 
Both  the  velocity  and  heading  channels  initially  display  unacceptable  transient  oscillations 
in  their  time  response.  The  control  system  designed  with  the  previous  aircraft  models  no 
longer  produce  satisfactory  results  with  the  new  aircraft.  A  fine  tuning  of  the  controller 
gains  results  in  acceptable  performance  using  the  new  second-order  models. 
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8. 7  Conclusions 


The  following  conclusions  are  drawn  from  the  results  of  this  study. 

1.  An  analytical  analysis  of  the  formation  flight  control  problem  revealed  the  necessary 
and  sufficient  condition  of  integral  feedback  in  order  to  control  a  formation  of  aircraft 
with  zero  steady  state  error.  Important  stability  conditions,  linking  the  aircraft  and 
formation  geometry  parameters  were  derived. 

2.  A  linearized  MIMO  plant  model  mathematically  describing  the  formation  flight  con¬ 
trol  problem  can  be  used,  through  pole  placement  techniques,  as  an  valuable  tool  in 
selecting  the  optimal  PI  controller  gains.  The  presence  of  non-linear  elements  in  the 
aircraft  models  and  kinematic  equations  requires  the  use  of  a  trial  and  error  method 
to  fine  tune  the  controller  gains. 

3.  The  formation  flight  control  system  developed  using  the  first  order  aircraft  models, 
can  successfully  track  formation  maneuver  commands  using  the  new  second-order 
models,  after  adjustment  of  both  the  PI  controller  and  linear  mixer  gains.  The 
slower  second-order  models  exhibit  a  more  realistic  time  response  by  not  immediately 
saturating  the  non-linear  rate  limiters. 

4.  The  formation  flight  control  system  designed  to  control  single  wing  aircraft  forma¬ 
tions  can  be  successfully  used  to  control  multiple  wing  aircraft  formations.  As  long 
as  the  formation  separation  distances  are  established  in  accordance  with  the  tracking 
error  of  the  lead  aircraft,  no  communication/coordination  between  the  wing  aircraft 
is  needed. 

5.  The  display  of  the  lead  aircraft’s  trajectory  in  the  wing  aircraft’s  rotating  reference 
frame  provides  a  graphical  means  of  evaluating  the  performance  of  the  formation 
flight  control  system.  Similar  to  Lissajous  figures  used  to  display  two  dimensional 
oscillation,  this  trajectory  is  used  to  evaluate  the  tightness,  speed  and  tracking  error 
of  the  control  system. 
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8.8  Recommendations  for  Further  Study 

Additional  research  in  the  area  of  formation  flight  control  is  needed  in  several  areas. 
This  work  will  in  part  address  the  limitations  and  assumptions  made  for  this  research. 
These  areas  are  described  below. 

1.  An  analytical  analysis  of  the  formation  control  design  problem  using  second-order 
models  should  be  accomplished. 

2.  More  extensive  computer  simulations  using  the  second-order  models  should  be  ac¬ 
complished  in  order  to  further  validate  the  formation  control  system. 

3.  Additional  research  should  be  done  on  the  multiple  wing  formation  flight  control 
problem. 

4.  A  method  of  quantifying  the  effect  of  the  non-linear  rate  limiters  on  the  control 
system  response  in  order  to  eliminate  the  trial  and  error  proceedure  used  to  select 
the  controller  gains  should  be  accomplished. 

5.  Actual  sensor  models  should  be  used  in  the  formation  control  system  to  allow  evalua¬ 
tion  of  particular  sensors  in  particular,  the  measurements  update  rate  is  an  important 
parameter. 

6.  Noise  should  be  incorporated  into  the  sensor  model  to  evaluate  its  interaction  with 
the  controller’s  gains  and  its  effect  on  formation  stability. 

8.9  Summary 

The  results  of  this  research  show  that  a  PI  controller  developed  for  automated  maneu¬ 
vering  flight  can  successfully  maintain  a  multiple  wing  aircraft  formation  without  collisions 
among  the  aircraft.  An  analytical  analysis  of  the  formation  control  design  problem  affirmed 
the  requirement  of  integral  feedback  to  control  a  formation  of  aircraft  with  zero  state  error. 
A  control  system  can  be  designed  to  conserve  the  energy  of  the  wing  aircraft  during  large 
heading  changes.  Finally,  second-order  transfer  functions  better  model  the  C-130  aircraft. 
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Appendix  A.  Second-Order  Aircraft  Models  and  Large  Aircraft  Flying 

Qualities 

An  integral  part  of  the  formation  flight  control  problem  entails  the  adequate  modeling 
of  the  formation  aircraft  dynamics  (see  also  (10)).  These  aircraft  models  are  based  on  wind 
tunnel  data,  actual  flight  test  data,  or  high  fidelity  computer  simulation  of  the  aircraft 
equations  of  motion.  For  the  purpose  of  formation  flight  control,  a  balanced,  accurate, 
yet  simplified  representation  of  each  formation  aircraft  is  needed.  Thus,  a  mathematical 
method  of  modeling  the  aircraft  dynamics  as  a  low  order  transfer  function  is  desired. 

The  technical  report  titled  “Second  -  Order  System  Models  of  High  -  Order  Plants,” 
by  M.  Pachter  J.  J.  D’Azzo  and  L.  E.  Buzogany  presents  an  analytical  approach  in  modeling 
high  order  dynamical  systems  (6).  The  response  of  the  system  to  a  unit  step  function  input 
is  used  as  a  means  of  evaluation.  The  system,  in  this  case,  is  the  actual  aircraft  control 
system,  subjected  to  either  a  step  velocity,  heading,  or  altitude  change. 

The  order  of  the  transfer  function  model  is  determined  by  the  dynamic  behavior  of 
the  aircraft  response.  Large  aircraft  exhibit  a  time  lag,  or  physical  delay,  in  response  to 
a  step  input.  A  first-order  transfer  function  is  inadequate  in  modelling  the  time  delay, 
since  it  includes  only  an  initial  positive  rate  of  change.  A  second-order  transfer  function 
is  desirable,  since  it  incorporates  both  a  finite  rate  of  change  and  a  “time  delay”  (6)[l]. 
The  next  two  sections  present  the  theoretical  background  and  procedure  used  to  derive 
equivalent  second-order  aircraft  models. 


A.l  Theoretical  Background 

First  and  second-order  system  responses  to  a  step  function  input  differ  in  one  signif¬ 
icant  aspect.  The  initial  response  of  the  first-order  system  has  a  positive  slope  at  f  =  0, 
while  the  second-order  system  has  invariably  a  zero  slope  at  t  =  0  (6)[l].  An  overdamped 
second-order  system  is  given  by. 


y{s)  _  ab 
u{s)  (s  -I-  a)(s  -i-  6) 


(A.l) 
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where  the  parameters  a,  6  €  R+.  Figure  A.l  displays  the  response  of  an  overdamped 
second-order  system  to  a  unit  step  input. 
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Figure  A.l.  Overdamped  Second  Order  Step  Response,  a=2,  b=l 


The  quantities  labeled  Onset  Delay(OD),  and  Onset  Rate(OR),  completely  determine 
the  response  of  the  second-order  system.  The  Onset  Rate  is  the  majcimum  slope  of  the 
system  response.  The  Onset  Delay  is  determined  by  the  point  of  intersection  of  the  time 
axis,  and  a  line  tangent  to  the  curve  at  the  maximum  rate.  The  values  of  OR  and  OD, 
can  be  empirically  determined  from  flight  performance  data,  or  they  can  be  analytically 
determined  using  analysis  (6)[7].  The  product  of  OR  and  OD  of  second-order  systems  of 
the  form 

s2-|-2Cu;„s-|-a;2 
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where 


Wn  =  undamped  natural  frequency 
^  =  damping  ratio 

results  in  a  universal  function  f{x),  as  defined  on  the  next  page. 

f{x)  =  ODxOR  (A.3) 

=  x^(l+- - - — Inij-l,  0  <  I  <  1.  (A.4) 

\  X  1  —  a;  / 

where 

X  =  — .  (A.5) 

a 

Note  that 

/(O)  =  0  (A.6) 

/(I)  =  ^-1«0.1036  (A.7) 

The  function  f{x),  and  parameter  x,  are  helpful  in  determining  the  values  of  a  and 
b.  Given  values  of  OR  and  OD,  the  parameter  x  is  determined  from  a  solution  of  Eq  (A.4), 
or  graphically  using  Figure  A.2.  Once  x  is  determined,  Eqs  (A.8)  and  (A.5)  are  used  to 
determine  a  and  6. 

a  =  ORx^  (A.8) 

An  important  restriction  on  the  determination  of  a  second-order  approximation  is  the 
maximum  value  of  /(x)  =  OR  x  OD,  as  defined  in  Eq  (A.7).  If  /(x)  exceeds  this  value, 
than  a  second-order  approximation  is  no  longer  possible.  Obviously,  in  this  case,  the  data 
came  from  a  higher  order  system.  A  similar  analysis  can  be  performed  for  an  underdamped 
second-order  system  (6)[5-6]. 

The  next  section  presents  a  design  example  in  deriving  a  second-order  model  for  the 
altitude  response  of  a  C-130  aircraft. 
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Figure  A.2.  Overdamped  Function  OR  x  OD 


A. 2  Design  Example 

This  section  illustrates  the  application  of  this  method  to  design  a  second-order  trans¬ 
fer  function  approximation  of  the  altitude  response  of  the  C-130  aircraft.  The  altitude 
response  is  derived  using  a  computer  simulation  of  a  six  degree  of  freedom  linear  state 
space  model,  equipped  with  an  Altitude  Hold  Autopilot.  This  autopilot  is  described  in 
detail  in  Appendix  B. 

The  first  step  in  the  design  is  to  obtain  the  step  response  (see  Figure  A.3)  and  the 
determination  of  the  two  parameters,  OD  and  OR.  Since  an  analytical  computer  model  is 
available,  OD  and  OR  can  be  analytically  determined.  Thus,  OR  and  OD  respectively,  are 
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Table  A.l.  C-130  Altitude  Response  2nd  Order  Parameters 


Parameter 

Value 

f 

.8  seconds 

Kn 

4.08  feet 

OR 

.1044  feet/second 

OD 

.409  seconds 

OR  x  OD 

.0427  feet 

X 

.08 

ah 

1.625 

bh 

.13 

are  shown  in  Eqs  (A. 9)  and  (A.IO). 


OR 

OD 


OD 


(A.9) 

(A.10) 


where 


OR  =  onset  rate 
OD  =  onset  delay 

h{t)  =  altitude  response  of  the  C-130  to  a  step  input 
t*  =  time  of  the  maximum  rate  of  change  of  h(t) 


The  unknowns,  t*  and  h(t*),  are  determined  directly  from  the  time  response  plots  of  h{t) 
and  The  time  response  plots  for  the  altitude  channel  of  the  C-130  are  displayed 

in  Figure  A.3.  The  unknown  parameters  extrapolated  from  Figure  A.3,  are  displayed  in 
Table  A.l.  OR  and  OD  are  calculated  using  Equations  (A.9)  and  (A.10). 


The  product,  OR  x  OD,  is  checked  against  Elqs  (A.6)  and  (A.7)  in  order  to  justify 
the  use  of  the  second-order  approximation.  The  value  of  x  =  ^,  is  determined  analytically 
using  Eq  (A.4)  or  graphically  using  Figure  A.2.  Elqs  (A.8)  and  (A.5)  are  used  to  determine 
Of,  and  bfi.  The  resulting  second-order  tramsfer  function  approximation  of  the  altitude 
response  of  the  C-130  is 


h{s)  .211 

hc{s)  ~  (s-|-.13)(«-l-  1.625) 


(A.ll) 
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Figure  A. 3.  C-130  Time  Response  Step  Altitude  Change 

The  actual  transfer  function  derived  from  the  linear  state  space  model  and  Altitude  Hold 
autopilot  is 


-.0007s‘2  -  .002s“  +  .025si°  +  .253s®  +  1.053s«  +  2.724s^+ 
h(s)  =  4.847s®  +  5,742s® +  4.725s‘*  +  2.19953+.596s2  +  . 081s +.004 

he(s)  .004si'‘  +  .051si3  +  .425si2  +  2.381s“  +  9.275s»®  +  25.85s®  +  52.883s«+ 

79.645s^  +  85.883s®  +  64.78s®  +  30.3s‘‘  +  8.77s®  +  1.456s®  +  .122s  +  .004 

(A.12) 

A  comparison  of  the  actual  altitude  response  and  the  second-order  approximation  is 
shown  in  Figure  A.4.  There  is  a  close  correlation  of  the  two  time  responses  with  maximum 
difference  of  about  10  feet  at  15  seconds.  The  second-order  transfer  function  does  an 
acceptable  job  in  modelling  the  actual  response,  while  providing  the  benefits  of  reduced 
order  modeling. 
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Figure  A.4.  C-130  Altitude  Time  Response  Compaxison 
A. 3  Summary 

The  use  of  accurate  yet  simple  models  of  the  formation  aircraft  is  a  key  factor  in 
the  solution  of  the  formation  flight  control  problem.  Second  order  transfer  functions  are 
desired  in  order  to  model  the  time  delay  present  in  most  physical  systems.  The  analysis 
presented  in  the  technical  report  “Second  -  Order  System  Models  of  High  -  Order  Plants,” 
by  M.  Pachter,  J.  J.  D’Azzo  and  L.  £.  Buzogany  contains  valuable  tools  in  the  development 
of  second-order  models(6).  These  tools  are  used  in  the  design  example  in  modelling  the 
14th  order  altitude  response  of  the  C-130  surcraft  as  a  second-order  transfer  function. 
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Appendix  B.  C-130  Aircraft  Models  and  Longitudinal  and  Lateral 

Autopilots 

An  important  factor  which  is  critical  to  the  design  of  a  formation  flight  control  system 
is  the  modelling  of  the  aircraft,  as  stated  in  Appendix  A.  The  aircraft  model  becomes 
integrated  into  the  system  plant,  and  thus  directly  impacts  the  response  of  the  system 
and  performance  of  the  flight  controller.  The  analytical  analysis  in  Chapter  VI  presents 
the  relationship  between  the  aircraft  parameters,  formation  geometry,  controller  gains  and 
stability.  Hence,  a  realistic  aircraft  model  is  required  to  have  confidence  in  the  results 
obtained  through  this  analysis  and  computer  simulations. 

A  second  consideration  in  developing  accurate  models  is  the  level  of  complexity 
needed  to  simulate  the  aircraft.  The  previous  research  used  first-order  models  based  on 
empirical  data  and  aircraft  time  responses.  These  models  represent  C-130  aircraft  equipped 
with  lateral  and  longitudinal  autopilots.  The  separate  velocity,  heading,  and  altitude  chan¬ 
nels  were  assumed  to  be  uncoupled.  While,  the  first-order  models  incorporated  an  initial 
rate  of  change,  they  could  not  simulate  the  onset  delay  present  in  many  physical  systems. 
Using  second-order  models,  both  the  onset  rate  and  onset  delay  are  incorporated.  Higher 
order  models  could  be  used,  but  the  increased  amount  of  analysis  and  computer  resources 
needed  would  oflfset  any  advantages  gained. 

Closed-loop  aircraft  models  incorporating  installed  autopilots  and  performance  data 
needed  to  derive  the  second-order  aircraft  models  are  difficult  to  acquire.  An  alternative 
approach,  for  use  in  this  study,  is  to  design  basic  autopilots  around  available,  bare  air¬ 
craft  state  space  models.  Computer  simulations  and  system  identification  techniques  are 
then  used  to  derive  the  reduced  order  models  used  in  the  formation  flight  control  system 
simulation.  The  latter  approach  is  used  to  develop  the  second  order  aircraft  models  for 
the  C-130  aircraft.  Linear  six  degree  of  freedom  state  space  aircraft  equations  supplied  by 
Lockheed  Aerospace  ace  used  (4). 

Basic  longitudinal  and  lateral  autopilots  are  designed  around  these  linear  mrcraft 
equations  using  root  locus  techniques.  The  time  responses  of  these  autopilots  are  then 
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analyzed  in  order  to  derive  decoupled  reduced  order  models  of  the  aircraft.  Table  B.l  lists 
the  flight  conditions  and  key  aircraft  parameters  used  in  the  model  development. 

Table  B.l.  Aircraft  Flight  Conditions  and  Parameters 


Parameter 

Value 

Description 

Type 

C-130H 

Aircraft  Type 

GWT 

12,000  lbs 

Gross  Weight 

ALT 

500  feet 

Altitude 

V 

306  fps 

Trim  Velocity 

The  complete  linear  six  degree  of  freedom  state  space  model  data  supplied  by  Lock¬ 
heed  is  attached  in  Appendix  A.2.  The  following  assumptions  are  used  to  simplify  the 
design  of  the  autopilots: 

1.  The  aircraft  body  axis  is  assumed  to  be  closely  aligned  with  the  stability  axis  for 
small  values  of  a,  angle  of  attack,  and  sideslip(l)[604-611]. 

2.  Velocity,  heading  and  altitude  step  inputs  are  separately  applied  to  the  closed-loop 
aircraft  in  order  to  identify  the  individual  transfer  function  models. 

3.  The  autopilots  are  optimized  for  command  inputs  of  10  fps,  10  degrees,  and  100  ft 
for  the  velocity,  heading  and  altitude  channels  respectively. 

In  the  following  sections,  the  design  of  the  C-130  autopilots  and  model  development  is 
presented. 

B.l  Longitudinal  Autopilot  and  Reduced  Order  Models 

A  Mach-Hold  autopilot  is  designed  as  shown  in  Figure  B.l.  Pitch-rate  feedback  is 
used  in  the  inner-loop,  while  proportional  plus  rate  feedback  is  used  for  the  outer  loop. 
The  gains  are  optimized  for  proper  damping  and  output  time  response.  The  autopilot  time 
responses  for  a  10  fps  input  is  shown  in  Figure  B.2. 

The  velocity  response  u  of  the  Mach-Hold  autopilot  closely  resembles  a  first  order 
response(l)[101].  In  this  case,  a  first-order  transfer  function  best  models  the  velocity 
channel  of  the  C-130  aircraft.  Figure  B.3  compares  the  first-order  approximation,  with  the 
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actual  autopilot  velocity  response.  These  plots  shows  that  the  first-order  transfer  function 
shown  in  Eq  (B.l),  very  accurately  models  the  mrcraft  performance. 
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Figure  B.3.  Time  Response  Comparison  of  Mach-Hold  Autopilot  and  First-Order 
Approximation 


«re/(s)  S  -1-  Tv  5  +  .1  ^ 

An  Altitude-Hold  autopilot  is  designed  as  shown  in  Figure  B.4.  The  aircraft’s  flight 
path  angle,  and  thus  altitude,  is  commanded  using  throttle  control.  A  second-order  cascade 
compensator  is  designed  to  stabilize  the  inner  loop  and  to  provide  sufficient  damping. 
Proportional  plus  rate  feedback  is  again  used  to  provide  the  best  altitude  tracking.  The 
Altitude- Hold  autopilot’s  response  to  a  100  feet  step  altitude  change  is  shown  in  Figure  B.5. 
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Figure  B.5.  Time  Respouse  of  the  Altitude-Hold  Autopilot  to  a  10  feet  Step  Altitude 
Command 


The  altitude  response  of  the  autopilot  is  clearly  second-order  or  greater.  An  over¬ 
damped,  second-order  approximation  of  the  altitude  response  is  developed  using  techniques 
described  in  Appendix  A.  Table  B.2  list  the  parameters  used  to  aetermine  the  second-order 
model. 


Table  B.2.  Altitude-Hold  2nd  Order  Parameters 


Parameter 

Value 

Onset  Rate 

10.436  feet /second 

Onset  Delay 

0.409  seconds 

Oft 

1.625 

bk 

0.13 

The  second-order  transfer  function  approximating  the  time  response  of  the  Altitude- 
Hold  autopilot  is 
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0.211 


Ohbh 


h{s)  _ 

^re/('S)  +  0’h){s  +  bh) 


(s  +  1.625)(s  +  0.13) 


(B.2) 


A  comparison  of  the  altitude  time  response  from  the  Altitude- Hold  autopilot  and  the 
second-order  approximation  is  shown  in  Figure  B.6.  This  figure  shows  a  close  correlation 
between  the  two  time  responses. 


Figure  B.6.  Time  Response  Comparison  of  Altitude-Hold  Autopilot  and  Second  Order 
Approximation 


B.2  Lateral  Autopilot  and  Reduced  Order  Model 

A  lateral  directional  control  autopilot  is  designed  as  shown  in  Figure  B.7(l)[176- 
188].  Yaw  rate  feedback  and  a  washout  filter  are  used  to  dampen  the  dutch  roll.  Sideslip 
feedback  is  used  to  provide  rudder  coordination.  Since  a  bank  angle  command  is  used  to 
turn  the  aircraft,  a  roll  angle  control  system  is  designed  around  the  coordinated  aircraft. 
A  limiter  is  used  to  limit  the  bank  angle  commanded  in  response  to  a  heading  error.  The 
outer  loop  utilizes  proportional  feedback  in  order  to  command  the  heading  error. 

The  response  of  the  directional  autopilot  to  a  10  degree  heading  change  is  shown  in 
Figure  B.8.  The  sideslip,  roll,  and  actuator  deflections  are  all  within  reasonable  limits. 
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Figure  B.7.  Lateral  Directional  Autopilot  Control  System  Design 


There  is  an  obvious  delay  in  the  heading  response  of  the  aircraft  due  to  the  necessity  of 
rolling  a  relatively  large  aircraft. 

The  heading  response  of  the  Lateral  Directional  autopilot  clearly  exhibits  character¬ 
istics  of  a  high  order  system.  A  second-order  transfer  function  is  designed  to  best  fit  this 
response.  While  there  is  not  an  exact  correlation  of  the  two  responses,  the  second-order 
transfer  function  better  models  the  transient  characteristics  than  a  first-order  transfer  func¬ 
tion.  The  second-order  model  is  shown  in  Eq  (B.3),  and  the  heading  response  comparison 
in  Figure  B.9. 


_  0.296 

V’re/(s)  ~  (^  +  +  b^)~  {s  +  0.544)(s  +  0.544) 
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Figure  B.8.  Time  Response  of  the  Lateral  Directional  Autopilot  to  a  10  Degree  Heading 
Command 


Figure  B.9.  Time  Response  Comparison  of  Lateral  Directional  Autopilot  and  Second- 
Order  Model  Heading  Response 
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B.3  Summary 


Lateral  and  longitudinal  autopilots  are  designed  around  linear,  six  degree  of  freedom 
state  space  models  of  the  C-130  Aircraft.  The  time  responses  of  these  autopilots  to  step 
command  inputs  are  used  to  derive  decoupled,  reduced  order  models  of  the  aircraft.  These 
models  are  used  in  the  continued  development  of  a  formation  flight  control  system.  Equa¬ 
tions  (B.4)  through  (B.6)  summarize  the  final  transfer  functions  used  to  model  the  C-130’s 
velocity,  altitude  and  heading  channels,  respectively.  Fiqure  B.IO  contains  the  final  system 
build  diagram  including  the  three  autopilots  for  reference. 


V{s)  0.1 

V,is)  ~  (s  +  0.1) 

h{s)  0.211 

h,(s)  ~  (s-f- 1.625)(5-f  0.13) 

V)(s)  _  _ 0.296 _ 

^e(^)  (5  +  0.544)(s  -b  0.544) 


(B.4) 

(B.5) 

(B6) 
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Figure  B.IO.  MatrixX  System  Build  Diagram  of  C-130  Autopilots 


B-10 


Bibliography 


1.  Blakelock,  John  H.  Automatic  Control  of  Aircraft  and  Missiles  (Second  Edition). 
John  Wiley  and  Sons,  Inc.,  1991. 

2.  Dargan,  John  L.  Proportional  Plus  Integral  Control  of  Aircraft  for  Automated  Ma¬ 
neuvering  Formation  Flight.  MS  thesis,  Air  Force  Institute  of  Technology,  Wright- 
Patterson  AFB  OH,  December  1991. 

3.  D’Azzo,  John  J.  and  Constantine  Houpis.  Linear  Control  System  Analysis  and  Design. 
McGraw-Hill,NY  1988. 

4.  McCarty,  J.M.  and  F.B.  Green,  “C-130,  LINEAR  SIX  DEGREE  OF  FREEDOM 
STATE  SPACE  MODELS.”  LOCKHEED  AERONAUTICAL  SYSTEMS  COM¬ 
PANY,  Marietta,  GA  1992. 

5.  Meyer,  R.T.  and  E.G.  Husband.  “Large  Aircraft  Flying  Qualities  Revisited,”  AIAA 
paper,  (90-2847):456-463  (1990). 

6.  Pachter,  Meir  J.,  et  al.  “Second  -  Order  System  Models  of  High  -  Order  Plants.” 
Submitted  for  Publication,  September  1992. 

7.  Porter,  B.  and  A.  Bradshaw.  “Design  of  Linear  Multivariable  Continuous- Time  Track¬ 
ing  Systems,”  Int.  J.  Syst.  Sci,  5,  1155-1164  (1974). 

8.  Porter,  B.  and  A.  Bradshaw.  “Singular  Pertubation  Methods  in  the  Design  of  Track¬ 
ing  Systems  Incorporating  Inner-Loop  Compensators  and  High-Gain  Error- Actuated 
Controllers,”  Int.  J.  Syst.  Sci,  12,  1193-1205  (1981). 

9.  Rohs,  Paiil.  A  Fully  Coupled,  Automated  Formation  Control  System  for  Dissimilar 
Aircraft  in  Maneuvering,  Formation  Flight.  MS  thesis.  Air  Force  Institute  of  Tech¬ 
nology,  Wright-Patterson  AFB  OH,  March  1991. 

10.  Trosen,  Dennis.  Aerial  Refueling  Control  System.  MS  thesis.  Air  Force  Institute  of 
Technology,  Wright-Patterson  AFB  OH,  December  1992.  In  preparation. 


BIB-1 


Vita 


Captain  Louis  E.  Buzogany  was  born  in  Cleveland,  Ohio  and  attended  Lake  Catholic 
High  School  in  Menter,  Ohio.  He  attended  Ohio  State  University  and  recieved  a  Bachelor 
of  Science  Degree  in  Electrical  Engineering.  Upon  graduation,  he  was  commissioned  a 
Second  Lieutenant  in  the  United  States  Air  Force  and  was  assigned  to  the  1000th  Satellite 
Operations  Group,  Offutt  AFB,  NE  from  May  1986  to  May  1991. 

While  at  Offutt  AFB,  he  performed  duties  as  a  Satellite  Mission  Planner  and  was 
responsible  for  planning  satellite  missions  for  the  Defense  Metereological  Satellite  Program. 
In  May  1991  he  was  selected  to  attend  the  Air  Force  Institute  of  Technology,  Wright 
Patterson  AFB,  OH  in  persuance  of  a  Master  of  Science  Degree  in  Electrical  Engineering. 


Permanent  address:  7101  Sandpiper  Ct. 

Painesville,  Ohio  44077 


VITA-1 


REPORT  DOCUMENTATION 

PAGE 

j  1.  AGENCY  USE  (Eo.iv.'  birnk) 

i 

2  REPORT  DATE 

15  December  1992 

3 

REPORT  TYPE  AN  ;■  OAT'S  TT  ViMrO 

Master’s  Thesis 

i  a.  Tiile  an'J  subtitle 

5 .  r  U  t--  ■  !>  S B  -  P  S 

AUTOMATED  CONTROL  OF  AIRCRAFT  IN  FORMATION  FUGHT 


6.  AUThOR(S) 

Louis  K  Buzogany 
Captain.  USAF 

f.  di?''~NT;7rnoM  and  addss;SS(F5; 

Air  Fence  Institute  of  Technology 
WPAFB  OH  45433-6583 


AF1T/GE/ENG/92I>07 


9  SPONSORIi9t  AGtNCY  NAMc(S}  AND  ADDrtFSSitS)  G;,  NO  MONilOiGNG 

A i~  f  \C  ''  K  r  r*’ ^  S: ;  ^  Vi ?■  f  « 

Wrigiht  Laboratories. 

HGS. 

Wrigbt-Pattetson  AFB.  OH  45433 


11.  SUPP*tF‘VliNTARY  nT)TcS 


")!«7'P!>i|Cio™’'AyA'LABILIYY  STATEMtNT 

Distribution  unlimited 


In  this  study,  an  automated  fonnation  cmtrol  system  for  an  aircraft  fonnation  comprised  of  one  lead  and  multiple  wing  aircraft 
is  analyzed.  Second-order  models  of  the  C-130  aircraft  are  developed  in  order  to  accurately  model  the  flying  qualities  of  large 
aircraft.  This  automated  formatioo  control  system  is  capable  of  controlling  the  C-130  aircraft  in  maneuvering  formation  flight 
thus  reducing  the  wing’s  pilot  wokload.  During  fonnation  flight,  flie  wing  aircraft  contiimously  measures  die  aircraft’s 
relative  position  witii  an  ideal  on-board  positirm  sensor.  This  information,  in  addition  to  Proportional  Phis  hitegral  feedback 
crmtrol.  is  used  to  maintain  the  aircraft  in  formatirm.  The  control  of  each  wing  aircraft  is  assumed  to  be  inA»pwi<twit  nf  nriiw 
wing  aircraft.  Other  than  for  nominal  formation  sqiaration  commands,  no  continuous  communication  is  ««nnne<t  between  die 
framation  aircraft.  An  analytical  analysis  of  the  fonnation  control  problem  reveals  diat  integral  control  is  needed  to  achieve 
zero  steady  state  error  in  the  sqiaration  distances  (after  a  formation  maneuver  is  executed).  This  conclusion  is  confirmed  using 
computer  simulations.  An  analytical  mediod  of  selecting  die  Proportional  Phis  hitegral  parameters  is  developed  by  identifying 
the  dominant  system  dynamics  and  residues  of  die  stq>  response.  In  an  attempt  to  reduce  die  fuel  consunqition  of  die  wing 
aircraft  during  formation  heading  change  maneuvers,  an  alternate  control  system  is  designed  to  conserve  the  energy  cS  die  wing 
aircraft.  The  resulting  automated  fnmiarinn  enntml  system  ftfifiectivriy  maintaiins  Ihf.  finrmatinn  tif  ahmift  thmiigh  a  <;rninhinf|tion 
of  velocity,  heading,  and  altitude  changes.  There  is  zero  steady  state  error  for  all  maneuver  and  sqiaratian  Higtanra  change 
responses.  Additionally,  the  sqiaration  transients  are  such  that  no  collisians  occur  between  die  aircraft.  Both  analytical  eolations 
_ and  computer  simulation  generated  time  responses  are  provided  for  reference. 
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